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PREFACE 


This manual contains solutions to the problems in Fundamentals of Statistical and Thermal 
Physics , by F. Reif. The problems have been solved using only the ideas explicitly presented 
in this text and in the way a student encountering this material for the first time would 
probably approach them. Certain topics which have implications far beyond those called for in 
the statement of the problems are not developed further here. The reader can refer to the 
numerous treatments of these subjects. Except when new symbols are defined, the notation 
conforms to that of the text. 

It is a pleasure to thank Dr. Reif for the help and encouragement he freely gave as this 
work was progressing, but he has not read all of this material and is in no way responsible for 
its shortcomings. Sincere thanks are also due to Miss Beverly West for patiently typing the 
entire manuscript. 

I would greatly appreciate your cal ling errors or omissions to my attention. 


R. F. Knacke 



CHAPTER 1 


1.1 


Introduction to Statistical Methods 


There are 6*6*6 = 216 ways to roll three dice. The throws giving a sum less than or equal to 6 


Throw 

1,1,1 

1,1,2 

1A,3 

1,1,2* 

1,2,2 

1,2,3 

2,2,2 

No. of 
Permutations 

1 

3 

3 

3 

3 

6 

1 


Since there are a total of 20 permutations, the probability is ^ 

1.2 


(a) Probability of obtaining one ace = (probability of an ace for one of the dice) x (probabil- 
ity that the other dice do not show an ace) X (number of permutations) = 

= (|) 5 = . 1*02 . 

(b) The probability of obtaining at least one ace is one minus the probability of obtaining 

none, or ^ 

1 - (§) = *667 

(c) By the same reasoning as in (a) we have 

1 2 £- 4 g, 

(£) (f) 57ir = - 040 



1*3 


The probability of a particular sequence of digits such that five are greater than 5 and. five are 

1 5 i 5 

less than 5 is (-) (— ) . Then multiplying by the number of pe imitations gives the probability 


v 2 y '2 
irrespective of order 


(W(h' - 2 46 
5151 '2' ~ 


1.2*. 


(a) To return to the origin, the drunk, must take the same number of steps to the left as to the 
right. Thus the probability is 





where N is even. 

(b) The drunk cannot return to the lamp post in an odd number of steps. 


1 




1.5 


r « 

(a) (|) 

(b) (Probability of being shot on the I-J trial) = (probability of surviving H-l trials) X 

th 5 N_1 1 

(probability of shooting oneself on the N trial) = (g-) (g-) . 

(c) 6 


-v ~ S r N 

m = nr = 0 since these are odd moments. Using n = (p (p+<l) 

and m 2 = (2n - u) 2 = U n 2 - 4nU + N 2 

i? = (2n - N) 4 = l6 i? - 32 Nn 3 + 2k N 2 n 2 - 8N 3 n + N* 
ve find m 2 = N and i? = 3U 2 - 2N. 


The probability of n successes out of N trials is given by the sum 

2 2 2 

W(n) = E E ... E w.w. ... w 
i=l j=l m=l 1 J m 

with the restriction that the sum is taken only over terms involving w^ n times. 

2 2 2 

Then V(n) = E w. E w, ... E v . If we sum over all i ... m, each sum contributes 
• i=l 1 j=l J m=l m ~ 

(wi+v 2 ) and 

W'(n) = (w x + w 2 ) N 

U i 

V ' (n) = E 7 /v* — r-r w., n w.^ n by the binomial theorem . 

^ 1 n=0 n:(w-n): 1 2 

Applying the restriction that w^ must occur n times we have 

> NJ n „ N-n 

= nl(N-n)i ' V 1 V 2 


Ve consider the relative motion of the two drunks. With each simultaneous step, they have a 
Probability of l/b of decreasing their separation, l/k of increasing it, and l/2 of maintaining 
it by ta kin g steps in the same direction. Let the number of times each case occurs be n^, n 5 , 
and n , respectively. Then the probability of a particular combination in N steps is 

t n l n 2 n 3 

B(a i n a n 3> ■ V n:ln 3 ; (|) ( 5 ) (|) W°3 - 11 


2 



The drunks meet if n^ = n^. Then the probability that they meet after N steps irrespective of 
the number of steps n^ which leave their separation unchanged is 

N n l n 2 


P = S „ n::nfc' ^ ^ 


n 3 =0 12 3 


where we have inserted a parameter x which cancels if n^n^. We then perfoim the unrestricted 

sum over n^ n^ n^ and choose the term in which x cancels. By the binomial expansion, 

1 1 i N , 2N / / 21T 

p ’ = x + + |) = (|) ( x ^ + x ) 


2U 2N 


Expansion yields 


® 2 hMtt <* 1/2 >° (^- 1/2 ) 


Since x must cancel we choose the terms where n = 2 U-n or n = H 
Thus P . (|) ( - K, l ' 

2 (mf 


(a) In (l-p) N_n = -p (U-n) ~ -Up n « U 


thus (l-p) N-n « e”^ 


00 InS}T = • • • ( N " n+ i) “ ^ 

r-\ Nl _n ,, _ % U-n U 11 n -Up X n 

(c) W w = niTu-nVi p (1 - p ) 15 ST P e - 51 


if n « If 


y \ n -X -X y \ n -X X 

S ^T e=e h-T=ee=l 


n=0 


» .n -X . v 

— - "A. o 

n=0 nI ^X 


y X -XX 

A — = X e e = X 
ni 


a 2 = 2 (B) 2 Z 4 - h 2 * X 

n=0 nI ^ n: 


/ \2 ? -o 

(An) = n - n~ = X 


(a) The mean number of misprints per page is 1 . 


.. n -1 

W(n) = — 7 e = — 7 

v ' ni ni 


W( 0 ) = e” i = .37 


3 



(*) 


2 -1 

P = 1 - E ^-r=.08 
n=0 n * 


1.12 


(a) Dividing the time interval t into small intervals At we have again the "binomial distribution 
for n successes in N = t/At trials 


n: 


W(n > = n X( lJ- n ) ~T P ^ 

X* X 
’ nJ 

tions in the interval of time. 


N-n 


In the limit At -> 0 , W(n) — e as in problem 1.9 where X = S f the mean number of disintegra- 

n« "t 


(b) W(n) = e' 4 


n 

0 

1 

2 

3 

4 

5 

6 

7 

8 

W(n) 

.019 

.076 

.148 

.203 

.203 

.158 

.105 

.061 

.003 


1.13 

2 2 

We divide the plate into areas of size b . Since b is much less than the area of the plate, 
the probability of an atom hitting a particular element is much less than one. Clearly n « N 
so we may use the Poisson distribution 


W(n) - fj e- 6 


n 

o 1 

3 

— r 

Hil 

.003 1 

.086 

.162 


1.14 


We use the Gaussian approximation to the binomial distribution. 

v 2 

1 


W(n) = 


(215-200)* 

2(400/1+) 




= .013 


1.15 


The probability that a line is in use at any instant is l/ 30 . We want II lines such that the 
probability that N+l or more lines are occupied is less than .01 when there are 2000 trials 
during the hour. 


IT 




.oi = i - Z 


2a 


n=0‘/^r a 


where 


- 2000 ~ ^ 
n = = 66.67, 


a = 




15 


= 8.02 





The sum may he approximated by an integral 


.99 = 




-§■ %/&? a 


exp [-(n-n) 2 /2a 2 ] dn 


p(l 2y 2 

• 9 9 = 1 - 7 =^ e ' dy after a change of integration variable. 


4 ~ M / V2 7T 

(-n-gj/CT 

The lower limit may he replaced hy -°o with negligible error and the integral found frcm the 
tables of error functions. 

We find * !Z n l i = 2.33 


N = 66.67 + (2.33) (8.02) - .5 * 85 


1.16 


(a) The probability for N molecules in V is given by the binomial distribution 


N I N v N_N o 

= NI(N -N)i ^ 

K o ' o o 


Thus 




0>) 


(N-fl ) 2 _ IT 2 - N 2 


o o 


1 - 


jj 2 


H 2 


H 2 


(c) If 


V « V Q , 


(AMP .. 1 
M 2 M 


(d) If V -» V Q , 


—9 

IT 


1.17 


Since 0 « — « 1 and since M q is large, we may use a Gaussian distribution, 
o 


P(tf)dN = 1 — i - exp [-(N-N) 2 /2 AN 2 ] dM 

(2?r AM 2 ) 


5 



1.18 


U N 

R = Z r. = -t Z §. 

_ i -a — 


— 1 


(S.’s are unit vectors) 


p2 ^ „ 2 . p2 v v* A „.p2 , p2 


r-r = tr Z s. + -t Z Z §.• §. =#+rZ Z cos e 

i 1 i / o ' X " J i / 0 

2 2 

The second term is 0 "because the directions are random. Hence R = Eft . 


1.19 


The total voltage is V = Z v. ; hence the mean square is 

i 1 


-s N — =■ N N _ 

v = Z v. + Z Z v. v. 

1 . / . X J 

1 1 / 0 


- 2 2 
where v. = vp and v. = v p. 
i x 


Hence 

1.20 


P = V 2 /r = (u 2 ^/ R) p 2 ri+(l-p)/lTp] 


(a) The antennas add in phase so that the total amplitude is E = ME, and since the intensity 

t 

is proportional to E^ 2 , 1^ = H 2 !. 

(h) To find the mean intensity, we must calculate the mean square amplitude. Since amplitudes 
may he added vectorially 

E , 2 = E, *E, = E 2 Z S 2 + E 2 Z Z S. *S, 
i i/O 

where the S^'s are -unit vectors. The phases are random so the second term is 0 and 


E 2 = HE 2 


or I t = HI 


1.21 


— g N ^ N N _ 

A = Z a + Z Za a 
n . n. . / . n. n. 
i i i / 0 i 0 


Since a = 0, the second term on the right is 0 and 
i 

i l 

.2 1 0.2 1 


(A 2 ) = H 2 (a 2 ) = H 2 10 3 a 


"but 


i_ 

(A 2 ) 2 = A„ = Ha c 


Equating these expressions we find N = 10°. 


6 



1.22 


N _ _ N 

(a) x =Es., but since = 'L, x = E -L = nL 


II N 


(x-x) = E (s i — ) + E E (s.--t)(s,- £) 

i i / j 1 0 


— 0**0 o 

Since (s -£) = £ - -t = 0, the second term is 0 and (x-x)^ = S a = Na . 


N _ II 

(a) The mean step length is -t. x = E s^ = E^= 

i i 


IT u ii 

(b) (x-x ) 2 = E (s. -£) 2 + E E (s.-l)(s.-l) , (s.~V)=-t--t = C> 

i j/i J 1 

To find the dispersion (s^-E ) , we note that the probability that the step length is between 

s and s + ds in the range -t-btO'L + b is • 


Hence 


7—^2 r Ul ° 'v 1 ) 2 ds i * 2 

- 4 , — ®— ‘ T 


, — \ 2 2 b 2 Nb 2 

(X-X) = E -r- = — 

i 3 3 


\1(9) d9 = H 


w(9) de = 


2 

27T a sin © d9 sin 9 d9 


(a) We find the probability that the proton is in 9 and 9 + d9 and thus the probability for the 


resulting field. 
From problem 1.2! 


W(9) d9 = 


sin 9 d9 


u 2_ ,\ , db 6u cos 0 sin 9 

Since b = >— (3 cos 0-1) we have — = -= 

a 3 a 3 


W(b)db = 2W(0)d0 = 


a 3 a/T db 


6a/ + |ia\ 


7 



(b) If the spin is anti -parallel to the field, 


W(b)db = 



Then if either orientation is possible, we add the probabilities and renormalize 


W(b)db = 



1.26 

0 “ . a« ° iks 

,W ■{ 48 ' S V ( s ) = TT dS 

Q may be evaluated by contour integration. For k > 0, the integral is evaluated on the path 



For k < 0 


e -kb 

residue (ib) = ^-r 


Q(k) = e -9 * 



8 



Thus Q(k) = e 

C/W -fj 


-kb 


,, -ikx -IT 
dk e Q 


'-Leo 


«-tr 


-ikx Hkb , 
dk e e + 


i. r 

ar Jr, 


dk e 


-ikx -Mkb 


-i r 

* Jn 


-Nkb , 
dk e cos k x 


This integral may he found from tables. 


(P(x) dx = - 


1 TTb dx 


7r x 2 + N 2 *) 2 

does not become Gaussian because the moments s n diverge. 

1.27 

From equation (1.10.6) and (1.10.8) we see that 

<?(*) = b X ** e_iix s. (k) -**V k) 


(i) 


where 


We expand (^(k). 


Qi(k) = 


ds e 


iks 


V-CO 


W.(s) 


Q.(k) = 1 + i s.k - | s* k 2 


In ^(k) = i Si k - | (As .) 2 k 2 


Q ± (k) = exp [i s.k - | (A Si ) 2 k 2 ] 

(P(x) = h: f dk exp[i(Z s ± -x)k - | Z (As i ) 2 k 2 ] 

“-<*> i i 


thus from (l) , 

This integral -may be evaluated by the methods of section 1.11. We find 

<P(x) = — exp [-(x-h) 2 /2ct 2 ] 

W2k a 


where 


N _ 

H=Zs. 


N 


c = Z (As . ) 
i 1 


1.28 


The probability that the total displacement lies in the range r to r + dr after IT steps is 

(p(v) d 3 r = J 2' f W(s x ) W(s 2 ) ... W(s^) d 3 ^ ... d 3 ^ 

( d3 l) 

where the integral is evaluated with the restriction 

IT 

r< Z s. < r + dr 
i=l 1 


9 



We remove the restriction on the limits by introducing a delta function 


CP (e) d 3 r = [ ... f W^) ... W(s^) S(r-Z s i )d 3 r 

O -CO 0 L 1 

ik* [Z s. -r] 

-j o ^ X ' 

In three dimensions 5(r - Z s ) = 5- / d^k e 

" i _1 (2tt) 3 -U " 


d 3 s. .. «d s„ 
-1 -N 


ik* [Z s^-r] 


thus 


or 


1.29 


<?(■*) = I'C'f W( -1 ) -*’ W( % ) ^3 X d3 - e 1 d3 ir-- d3 % 

Ptr) -^5 /_ i3 i e “ £l •••£ d V<%) 

<p (r) = — ^ d 3 k e — Q^(k) where Q(k) = d 3 s_ e"^*— W(s) 


i — *— N 


For displacements of uniform length hut of random direction W(s) = ( s ~ ) , in spherical 

% s 

coordinates. W(^) is a properly normalized probability density since 

f f f s 2 sin 0 ds d0 dcp = kn f = 1 

Jn Jn Jn W B 2 J n ^ 


1 0 J 0 u 0 

We find the Fourier transform Q(k) 

<Kk) = fff ^ e^w(s] = rrr 

J J J J Q J Qj Q 


> 27T p7T p“ 


s sin 0 ds d0 dcp 


S(s-t) £ iks cos 0 
% s 2 


r* sin 0 d0 ikL cos 0 sin kt 

mar J n —sr~ ’ -TT 


From 


we have, for N = 3 


(27T) 3 


.3, -ik*r -N/, \ 
d k e Q (k) 


1 2ir p7r p® 


(P(r) = — ~r / II k“ sin 0 dk d0 dcp e 

(27r) d ^0 


-ikr cos 0 sin 3 ki 

(kt) 3 


The angular integrations are trivial and 


?w ■ sk /> 


sin J kt sin kr 


To evalu a te the integral, we use the identities 


*3 

sin kE, sin kr 1 . . . . 0 ,1-cos 2 kEU 

5 = —= sin kr sin kc ( — — — 

k k 2 


10 




CHAPTER 2 



II 



2.3 


(a) Tiie probability of displacement x is the probability that cp assumes the required value in 
the expression x = A cos (wt + cp) . 

r(<p)dcp = and P(x)dx = 2v(cp) dtp = ^ 


w( 


7 r 


vhere the factor of 2 is introduced because two values of cp give the same x. 


Since 


we have 


dcp 

dx 


A sin (ut + cp) r . j p 

VA - x 

P(x)dx = 


dx 




(b) We ta kp the ratio of the occupied volume, 8v, in phase space when the oscillator is in the 
range x to x + dx and E to E + 8E, to the total accessible 
volume, v. 


The area, A g , of the ellipse is up^x^. 


Since 


2 2 2 
* _ P_. m 03 * 

E - a + 2 



we have 


Hence 



A e = 7r ™ ▼ = if 5E * 

' mw 


*(*)■* - a - 


V 27r o-r» 
~ 533 


For any x, p ^ = 8E 

, » , mu dx mu dx 

T - w j - g — 

V2mE - mux 

12 2 

where we have used (l). From (l) and x = A cos (ut + <p), it follows that E = — u mA . 


Thus 


2.k 


r(=0a* ■ i 7=§=r 

■UJ 


(a) The number of ways that N spins can be arranged such that n^ are parallel and n g are anti- 
parallel to the field, is 

HI 


r 'V ■ 


and 


“W - f (“i> ffe 


12 




Since 


Thus 


= -(n^) nH = -(a^-H) nH, ve have ^ = | (N - ^), n 2 = | (N + ^) 

„/ E \ I* M_ 

2 hE 

v 2 2nH y * '2 2pH' * 


(h) Using Stirling's approximation. In n! = n In n-n, ve have 

Ni 


In 


^T#3^yi ■ - n l ^ T + ( V B) ^ {1 ' T* 


a(E) - - it* - fe) ^ - & - §<» ♦ ♦ *> ♦ * & 


(c) We expand In f (n^) = In , (uln^) 1 a ' t>out tlie maximum, n^ 

In f(n 1 ) = In f^) + § B 2 


( 1 ) 


B 2 = 


In f(n 1 ) 


where ti = el. - n 

i. 1 ^ I , c. 

i “i 

- 1 |b 2 | t] 2 

From (1) f^) = f^) e 2 1 21 (2) 

f (n^) is evaluated by noticing that the integral of f(n^) over all n^ must equal the total number 


of permutations of spin directions , 2^. 

f 


‘I B 2 I" 


dij = 2 11 


*<=l) * 


n^ is found frcm the maximum condition 


d In f(n^) 


dn. 


= 0 = -In n^ + In (N-r^) , or = | 


and 


N - 

d 2 In f(n^) 


_ 1_ _ 1 

dni 2 

IV 

IV »v 

“l N -°l 


“l 


k 

N 


Substituting these expressions into (2), ve have 


\ -B /2 

2 
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Since 


“M * r <“i> §s > 


and 


n i ■ I (s - s ); 


it follows that 


a(E) = 


2.5 

(a) 

If is an exact differential. 

Since second derivatives are equal, 

(b) If f is exact, f 


,N 


y-rr g 


e *P f" 4 (of^) 3 


N v 2pE / J 2ioH 


^f = Adx + Bdy 

A.j£ B.g 

ox dy 

cA _ SB 
ciy Sc 

^ f = f (B) - f (A) =0 if A = B 


2.6 

(a) not exact 


d 2 f 


d 2 f 


O X / O X 

ST 6y ' dy ox 


(b) J dy + J (x 2 -2) dx = i, J" (x 2 -l)dx + J 2dx = ^ , J (x 2 -x)dx + J' x dx = 

(c) exact 


(a) 


.2,2 


0. f = 


1,1 


x + 


(2,2) 
xl (i.i) 


e" 

s 


= 1 


for all paths 


2.7 


(a) The particle in a state with energy E does work, li W, when the length of the box is changed 
to L + dL . This work is done at the expense of the energy, i.e., c^W = -dE. Since c^W = F dL , 

XX XX 


dE 


it follows that F = - — . 

x dL x 


x 2 ft 2 ✓ 

(b) The energy levels are given by E = ■ ■ ■ ( 


2 2 
n n n 

JE_ + _2_ + 

.2 2 + 2 

L L L 

x y z 

2J2 


- * 1 5E _ ^7t 2 h 2 / JC_ 

_ P ’ W-z * ' V* K ' "VA v L 2 

r \ X 


— %K g 
p mV 


s n 

(f 


where 


V = L L L 
x y z 


If L. . = 


x 

2 :2 


, = L = l , by symmetry n = n = n , and using the expression for E we have 

j ** x y z 


u> M 



2.8 




p = a V , a may "be found by evaluating this expression at a point on the curve. In units of 

6 2 o O 

10 dynes/m and lO^cm, a = 32. For the adiabatic process, 

AE = Eg - E a = -J 32 V -5 / 3 dV = -36 =-3600 Joules 
Since energy is a state function, this is the energy change for all paths. 

08 




dV = 22,400 joules 


Q = AE + W = -36OO + 22,400 = 11,800 joules 
P = 32^- 4jr (V-l) 

W =J 32 - 2 k (v-l)J = 11,550 joules 

Q = 11,550 - 3600 = 7,950 joules 

r.8 

W = J (1) dV = 700 joules 

Q = 700 - 3600 = -2900 joules 


15 



CHAPTER 3 


Statistical Thermodynamic s 


3-1 

(a) In eauilibrium, the densities on "both sides of the box are the same. Thus, in the larger 
volume, the mean number of Ne molecules is 750 and the mean number of He molecules is 75 • 


1000 , 100 


(1) 

3-2 

From problem 2.b- (b) 


p - (?) (?) - 10 


-185 


( a ) 


In fl(E) - - 2 (H - ^ In 2 (1 - " 2 (H + ^ In 2 ^ 1+ - ; ^ 

6E 


where we have neglected the small term In ^ - . 




1 - -(1 — — ) 
x 2 k H^H' 


Consequently E = -NpH tanh ^ 

(b) T < 0 if E > 0 

(c) M = n(n^-n 2 ) = n(2n^-Il) where n^ is the number of spins a.lllgned parallel to H. 

771 

— ) 


1 T7* 

Since ^ = i(N - from problem 2.k (a). 


M - u(H - jig - S) - - | 


M = Nu tanh 




it 


3-3 


(a) for system A, we have from 2.k (c) In 0 (e) = 


*2 


In 2 




2 t i 2 H 2 lT 


(1) 


ttN 

2 


where we have neglected the term In 


SE 

2uH 


In equilibrium. 


3 = §=r In fl(E) = 
3 = 3’ 


u 2 H 2 H * 


3’ = 


4’ 


2 H 2 H’ 


Hence 


u ,2 N* 


(b) Since energy is conserved, E + E* = bNpH + b’H’ioH 

Substituting from (a) we find E = + b'K'ji'H } 


(c) 


p 2 !! + h ,2 H* 

Q - E - Muz . MgO-VV* - V a ri free (2) 

HIT + li'TT 


( 2 ) 
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(d) P(E) = J1 (e) fl'(E 0 -E) where E q is the total energy, E^ = hWpTT + b'N'ii’H* 

E 2 


From (l) 


P(E) = exp 


2 2 
L 2^2 N 


exp 


1 

L , 2 H 2 tt' J 


L 2 s , 2 H 2 ir 

This expression may he revrritten in terns of the most probable energy E, equation (2). 

P(E)dS = C exp pCg rlL dE 

L 2 a J 


where 


? 2 2 

2 = u u ,g H g M t 




4 


The constant C is determined by the normalization requirement 

(E-E) 


L 


C exp 


2a 


dE = 1 


Thus 

(e) 

(f) 


3^ 


P(E)dE = -=y- exp 
a^T a 


(eJV 


2a 2 -1 


dE, where a is given above . 


(A*E) 


2 = 2 = u 2 u» 2 S 2 M , 
n 2 N+n , 2 N* 


From equation (3)> A*E = qE Vn" 

~ = 11 NfbuHW + TaVU'H) __ ,, Nu^g 


liN +11 


, 2 N’ 


hence 


A*S 

E 


4’ 


jiL 


N* » N 
N' » N 


n bi/iT 


(3) 


The entropy change of the heat reservoir is As' = " a^d for the entropy c hang e of the entire 

Q 

system, AS + AS* = AS - ^7 ^ 0. , by the second law. 


Hence 


AS a 


3.5 

(a) In section 2.5 it is shown that n(E)= V N x(e), where V is the volume and X (e) is independent 
of volume. Then for two non-interacting species with total energy E 

N l +N 2 


0 >) 


n(E) = C ^(E) fi 2 (E o -E) = CV X X 1 (E) X 2 (E) 

In n(E) = (M^+Ng) In V + In C + In X 1 (e) Xg(E) 
-IS 1 ( N l +N 2 ) 

P = f §v In 0(E) = | -- v — , or pv = (l^+N^kT 
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3-6 

2 atmospheres. 


CHAPTER b 

Macroscopic Parameters and their Measurement 

k.l 

(a) AS^ = me J ^ = 4 l 80 In = 1310 Joules A 

Since Sres = 

ASres = “ = - Tf (373-273) = -H20 joules A 

Aot = A + ^res - 190 JouleS A 

(h) The heat given off "by either reservoir is 5= mc(50) = 2.09 X 10^ joules. Since the 
entropy change of the water is the same as in (a) we have 

AS = 1310 - 2,09 = 102 joule s/°K 

(c) There will "be no entropy change if the system is brought to its final temperature by inter- 
action with a succession of heat reservoirs differing infinitesimally in temperature. 

k.2 

(a) The final temperature is found by relating the heats, Q, exchanged in the system. The heat 
required for melting the ice is m^-t where m^ is the mass of the ice. Since for the other compo- 
nents of the system, Q = mcAT, where m is mass and c is the specific heat, we have 


+ ^(melted ice) + ^(original water) + ^(calorimeter) ^ 

m i l + ^(V 273 ) + Vw(V 293) + m c C C ( T f‘ 293) = 0 

(30) (333) + (30)(^.l8)(T f -2T3)+(200)(4.l8)(T f -293)+(750)(.Ja8)(T f -293) =0 

T f = 283°K 
m.-t 

(b) The c hang e in entropy when the ice melts is For the other components 

r Tf dT T f 

AS = me / — = me In — 

J- i 

£ 

Hence AS=^ + m w cln||| + m i cln||| + m c c c In §|l = 1.6 JoulesA 
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(c) The work required is the amount of heat it would take to raise the water to this temperature. 

w = Vw <V*i> + m c c c 


( 750 ) (. 418 ) + ( 230 ) ( 4 . 18 ) 


(293-283) = 12,750 Joules 


4.3 


2 r Q = cdT+pdV by the first law 


= c dT + y dv 


Since pV = BT 


AS 


-r 


- 

T 


■T, 


dT ^ 

C T + 


r f Rdv 
Jv V 


T V 

f f 

= cln — + Rln— independent of process 
1 V i 


4.4 


The spins become randomly oriented in the limit of high temperatures. Since each spin has 2 
available states, for If spins we have 

S(T -»«) = k In n = k In 2^ = Nk In 2 (l) 

rp 

S(T)-S(0) = J dr 


But since 


we have 


S(T ~ m) '4* c i ( f ' 11 f ' c i (1 - ln 2) 


( 2 ) 


where we have set S(o) = 0 by the third law. From (l) and ( 2 ) Hk In 2 = c^(l-ln 2) 

Hk In 2 


C 1 ~ 1 - In 2 


= 2.27 Hk 


4.5 


The entropy at temperature T is found by evaluating S(t) - S(o) = f dT. Labeling the 

J n 


undiluted and diluted systems by u and d, we have 


S u< T l> - S u (0) =/^ c l( 2 - 1) f = c l( l-ln 2) 


S d (Ti) - S fl (0) = f 


**1 


T fl 

C 2 T 2 T = C 2 


S d (0) = S u (o) by the third law. Since there are only as many magnetic atoms in the diluted 
case as in the other, we must have S^T^) = .7 S^(T^) . Hence 


— = .7 (1-ln 2) = .214 
C 1 
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CHAPTER 5 


Simple Applications of Macroscopic Thermodynamics 


(a) pV 7 = constant for an adiabatic quasistatic process. Substituting pV = v KT ve have 

vR T.V 7 v ET.V 7 V_ 1-7 

- 1 — = f - f — or T = T ( — ) 

V i V f f i k ; 

(b) The entropy change of an ideal gas in a process taking it from T^V^ to was found in 

problem it.3 T V 

AS = vCy In — + vR In — 

AS = 0 for an adiabatic quasistatic process. Hence 

T V R//c V 
° = In^ ( ^ ) 

A i v i 

Since R/c„ = c -clVc^ = 7-1, and since In 1 = 0, it follows that 
* P Y V 1 --y 

V f 

T = T ( — ) 
f i k V ± ' 


(a) W = p d7 = area enclosed by curve = 31^ joules 

(b) Since the energy of an ideal gas is only a function of the temperature and since pV = KT, 

P V "P v 

A „ 3 „ , c c ‘A An 
AS = ^ AT = | R (— g-) 

= ^ (6-2) X 10^ ergs = 600 joules 

(c) Q = AE + W where AE is calculated in (b) since the energy change is independent of path. 


Again the work is the area under the curve. 


Q = 600 + 1+00 + 


1157 joules 


_ _ SE _ 5 R 

®V " cff " 2 R 


W = / p dV = area under curve = 1300 joules 


(c) To find the heat, Q, we first find the energy change in going from A to C. This, of 
course, is independent of the path. 

AE = c^AT = | R(T c -T a ) = | (P C V C - P A V A ) = 1500 joules 
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Then since Q = AE + W 


(d) From problem 4.3 


5.4 


Q = 1500 + 1300 = 2800 joules 


AS 


' C V 


In + R In 
A 


V 

c 



= f-R In 


12 X 10 2 /R 
6 X 10 2 /R 


+ R In 


3 x 10 3 
10 3 


= 2.84 R = 23.6 joules/°K 


(a) The total system, i.e. , water and gases, does not exchange energy with the environment, so 
its energy remains the same after the expansion. Since the water cannot do work, its energy 
change is given hy the heat it absorbs, AE = Q = C^AT, where is the heat capacity. The 
energy of the two ideal gases is, of course, dependent only on temperature, AE = GAT, where C is 
the heat capacity. 


“(vater) * “(He) + “(A) ' 0 


Sl< T rt> + C H e< T rt> + C A< T rt> ’ 0 


Thus T f = T^, no temperature change. 

(b) Let -t be the distance from the left and v the number of moles. 

V He P He V He = (5)(30) _ , 

V A " P A V A "tl K 50 T - 3 

Hence, in equilibrium -t = 3(80-t), -t = 60 cm. 

V f 

(c) Since there is no temperature change, the entropy change is vR In — . 


AS = A£> a + AS He = |ln|2 + Rln|^ = 3.24 joules/°K. 


5-5 

(a) The temperature decreases. The gas does work at the expense of its internal energy. 

(b) The entropy of the gas increases since the process is irreversible. 

(c) The system is isolated, Q = 0, and from the first law 

AE = -W = - (V f -V Q ) 

The internal energy of the gas is only a function of temperature , AE = vc^(t^,-T q ) 

Thus vCy(T f ^ o ) = - 2£ (V f -V Q ) (1) 

In equilibrium p = ^ and since pV^ = vRT^ , we have 
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VAST, 

V = = 

f mg 


1 r V 1 

Substitution in (l) yields T f = ± '+ ( R /^) [ T 0 + 


5.6 


The equation of motion is 
Since the process is approximately adiabatic 


mx = pA - mg - P q A 


and 


pV^ = constant = (p + 5£) v/ 
o A o 


AV ' 

** = (P Q + ?) ~ ^ - mg - P„A 


(Ax)' 


( 1 ) 


( 2 ) 


The displacements from the equilibrium position, — , are small so we can introduce the coor- 
V V A 

dinate change x = ^ + T) and expand about —■ . 




7 V 

X / O . \' o o 

(t + ’l) 


We keep only the first and second terns and substitute into (2) . 


mq 




This is the equation for harmonic motion and the frequency may be determined by inspection. 

1 

Hence , 


V = 2^ 


(P 0 + T?) 


A ' V m 
o - 


p p 

% v m V. 


P Q A + mg A 


5-7 


The volume element of atmosphere shown must be in 
equilibrium under the forces (pressure) x (area) and 
gravity. Then if n is the number of particles per un i t 
volume, m the mass per particle, and g the acceleration 
of gravity. 



z+dz 

z 


p(z+dz)A-p(z)A = -n(Adz)mg 

dz = - dz, where = m 
dz 3J, 7 w 


N, 
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Since p = nil. 


( 1 ) 



( 2 ) 

( 3 ) 



Consider the element shown. Ax is its initial width. The wave disturbance causes the left face 

to move a distance g and the right a distance |4Af. Let P q and p be the equilibrium pressure 

and density, and p = P-P o is the departure from the equilibrium pressure. 

1 AV _1_ AA|_ 

K S “ V P-P Q AAx p 


or 
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where we have neglected the pressure difference Ap across the slab, 
from the equation of motion. 


m 


->2 

— -§• = PA - (P-*AP) A = -APA 

dt £ 


Since 


we have 


Then 


hut since the 


AP 



Ax 

ap 


SE = 

Sc 

aA 2 

a? 

StJ 

Sc 


AAx 


= “ P 


A 

St 2 


« is. 

3 St 2 


SjN a 2 E _ 1 a 2 p 

Sx/ ^2 


at 


^x 2 


order of differentiation is Immaterial, 

a 2 p 1 a 2 p 

at 2 K S P ax 2 


5.9 


(a) For an adiabatic process pV^ = const. 

tt av v . i av 1 

dp yp S V dp 7P 

(t) u = (p k s )^ = ( 


The displacement is found 


Substituting 

p = ^ and p = ^ > ^e find 


u = 

(c) 

V * ) 

1 

u « T^ , independent of pressure 

(d) 

u = 35^ meters/ sec. 


5.10 


Substituting the given values into the equation 


C V = 


2 2 
VT a 


we find 


cl. = 2k. 6 joules mole “deg 1 , 


7 = l.lt 


5.11 

When the solid, of dimensions xyz, expands an increment dV, we have 


2k 



V+dV = (x+dx) (y+dy) (z+dz) 


To first order 


V+dV = V + yzdx + xzdy + xydz 


_ldV_ldx ldjr ldz 
“ " V dT x dT y dT z dT 


Therefore 


a = 


= 3 °L 


5.12 

For the adiabatic, quasistatic process 


33 - (M) ® + (l) T *=■ ■ ° 


Dstog (11 ‘f 


and the Maxwell relation to 




_/dv\ 
\Sv , 

\ / "D 


we have 




c 

aVT 


. 1 /oV 

where O' = ? ^ 


Thus 


AT = — Ad 

pc 


where c = -£• 

P pv 


5.13 

From the first law c^ = T 


Substituting the Maxwell relation 


l * T (i 

T v 

si ■ H/ 


fLX = T 






1^1 and the definition a = - — fvs=\ we have 


.S 




-T (^) ov = a 2 vT - vT g 


5,1k 

(a) 


TdS = dE - FdL 


(b) From (l) we may read off the Maxwell relation 


Ss 

5l 


oF 

Sr. 


Since 


f T C r 

(c) S (L o ,T) - s(L o ,T q ) = J -i 

o 


F = aT 2 (L-L o ) 

l) T - - 2aT HI 

_ t _ 

dT * = J ~ ^ dT* = b(T-T 0 ) 

T o 


(1) 

( 2 ) 

(3) 


25 



Hence 


S(L,T) - S(L q ,T) = J (^) - 2 aT(L*-L o )dL' = -aT(L-L 0 ) 2 


S(L,T) = S(L q ,T 0 ) + b(T-T 0 ) - aT(L-L Q ) 2 
(d) In tlie quasistatic adiabatic process = 0 

S(T o ,L o ) + t(T f -T o ) - aT ± (L f -L o ) 2 = S(T o ,L q ) + -aT f (L f -L Q ) 2 


Hence 


m _ m 

f i 


t-a(L ± -L q Y 
t-a(L f -L q Y 


Thus if L^> L a we have T f > i.e. the temperature increases, 
(e) From the first law />c\ 

^ ' T (s) t 



Substituting (2) and ( 3 ) we find 


vSirl = “ 2aT (L - L o } 


(?) S(L,T o )-S(L o ,T o ) 



p j. u_ ax * P 

S(L,T) - S(L,T 0 ) -/ -fr- -/ 


T C dT* r T hT* - aT'(L-L )' 


T* 


dT* 


Thus 

5.15 


(a) 

(h) From (l) 


= b(T-T o ) -a(L-L o ) 2 (T-T q ) 

S(L,T) = S(T o ,L q ) + b(T-T o ) - aT(L-L Q ) £ 


•d Q = TdS = dE - 2a £ dx 

as-f 

T T 


( 1 ) 


Hence 


dx + 


1 fo£\ 

dT = T 


, 1 

+ t VSF; 


_ 2ot , 
«3T - — dx 


Equating coefficients, we have 


: d _ - <1 


+ 2ffi 


26 



From (l), ve can read off the Maxwell relation 


SS\ /§(-2g£)\ do 


Since a = o q -off. 


(¥) 

W T 


= 2 lea + 2 la - 2ic£T = 2la 

W T 0 ° 


If the film is stretched at constant temperature 


E(x) - E(o) = 2 ^o q x 


(0 -> x) = - J' ~F dx = ~ J 2o& dx' = -2 <j£x 


By the first lav 


as -f + 

T T 


Hence 


SsA 3 . /SS N \ 1 /SE\ , . 1 /"dE\ . zfV 

H dv + W v ” ■ T <vV T 4V + T (Sf ) ® + — d ” 


Equating c oefficients, ve have 


II -*(g) 


From (l), ve may read off the Maxwell relation 


Ss'N . fbY' 

■ df {-m. 


= zfT 


If one mole is produced 


AE = zfT - zfV 


By equation (5-8.12) 


I) T ■ 1 (i) v - ’ 


Substituting p = n Iff [1 + B„(t) n] ve find 


= P + ^ -3S " P = z 2 *? 
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5.18 

(a) 


dE = 


i) T ® + (I) v 41 


@x-v-(i) =<i; 


from (5.8.12) 


(h) From the first lav, dE = TdS - pdV. Hence 


(c) For a van der Waals gas 


/ds\ „ /'oS\ p 

0 ' T (sf) E ' 5 “* (,Sv) E * t 


vRT v a 


dTN v a 


T ^ = r v 2 r&\ w = r v 2 dv = .avf a _ jj 

21 X V, We X V 2 °V ^ V 2 V l- 


(a) From (p + ^j(v-b) = RT 


& . SL -0 


T (v -b) 2 


XSSX c 6a 0 

W4“ (. -h)3 ' 7 1 * 

T v c ' c 

O T 

Solving, ve find a=jRtv.b=- 

o c c 3 

, RT 

(h) Substituting a and h in (l) yields p =4 — - . 

c 


P ’ + ^)( V ’ "t 


To find the inversion curve, ve must have 


'8T\ V /T roV\ ' 

,W H " c; (v <Sr> p - 1 




From the van der Waals equation 


V 

T 


dp 


RdT 

v-b 


(v-b) •<r > s 



thus 



On eliminating v and putting the equation in terns of the dimensionless variables of problem 


5.19# it follows that p 

p« = 9 _ 12 {spP' - V? ) 



5.21 



Substituting these expressions into ( l) we find 
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(b) Frcm (2) 


r T dT r* a'd-d* 

J T T = J 4 f C * 


o 1 + 


T = T exp 
o 


" 

f 15 a'dd 

” 

V. 

u , - 

o 1 + — 

r cT T 

V 


5.22 

The system may be represented by the diagram 

^2 


<l n 


W 


(a) By the second lav 
For TnBYiTmTTT' heat, the equality holds and since = q^“W , it follows that 


g l *2 
AS = ^ 


*1 T i 


O’) 

5-23 

We have the system 


W Tl -T o 

°i 

^ = 11.9 



(a) W = q 1 -q 2 = C^-Tj) + C(T 2 *T f ) = C^-KT,, - 2T f ) 

(b) 3y the second lav. 


( 1 ) 


Thus 


T f 

r 1 

/ 

J T 

i l 


AS § 0 

T 2 

so + r f s ®- clJ1 !s_ 
T 4, T T 1 T 2 


5 0 


T f 1 


It follows that 

(c) The maximum amount of work will be obtained when T^ = ^ T^T^ . 
(1) W = C^-fTg - 2T f ) = 
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To freeze an additional mass m of vater at T 0 , heat mL must he removed from the ice-water mixture 

resulting in an entropy change AS, = - jjp 1 . The heat rejected to the "body of heat capacity C 

f T f dT ^f 

increases its temperature to T„ with an entropy change AS = C L — = C In — . By the 

mL Tf u "o T o 

second law AS^ + ASg = - — + C In ^ S 0. For minimum temperature increase and thus minimum 

“ o o 

heat rejection the equality holds and it follows that 


Tf = T q e 


mL/T o C 


The heat rejected is Q = C(T^-T o ) = CT q (e 


mL/T C 
■ ° -1) 


The prohahility that the system occupies a state with energy E is proportional to the number of 
states at that energy, fi(s). By the general definition of entropy, S = k In J1 (e), we have for the 


prohahility 


P « n(E) = e 


_ sA 


To raise the weight a height x, an amount of energy, mgx, in the form of heat, Q, must he given 
off hy the vater. The temperature change is found from 

Q = CAT = -mgx 
AT = -mgx/c 

where C is the heat capacity of water. For a mole of liquid and x = 1 cm, AT is on the order of 
10 ^ °K. Hence we may approximate the water as a heat reservoir and find for the entropy when 


the weight is at height x 


5 = s o-r = s o- 

o 


where S q is the entropy at height x = 0. Thus (l) becomes 

-mgxAr„ 

P(x) = Ae 

s oA 

Here we have absorbed the term e into the constant A which is found to he kT /mg by the 
usual normalization condition. Then the probability that the weight is raised to a height L or 


more is 


r-\ w -TV i 

L pWax -/ T -£ 


oo hr -mgxAT 


dx = e 


-mgiAr 


Substituting L = 1 cm and the given values yields 


GW 10 
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5.26 

la the processes a -> h and c -> d no heat is absorbed, so by the first law W = - AE, and since 
a? = ry\T , where C is the heat capacity, ve have 

V* * -<V E a) * -*V 

"c^d ' -< E d- E =) * ' T d^e> 

■y— T 

However, in an adiabatic expansion TV' = const. 


Hence 


Vd, * ^ ( 1- iy - -S T b C 1 - <v^ ) 

w <~d - ^ T o - 1 ) * "V T = ( 1_ <^> ) 


The volume is constant in process b-»c so no work is performed. 

= (Eg-E^) = vGy. (T c ^) 


7-1 


W , + V j /i aX 

s — c — > cl _ f c\ 

11 = s. vv 

CHAPTER 6 

Basic Methods and Results of Statistical Mechanics 

6.1 

(a) The probability that the system is in a state with energy E is proportional to the Boltzman 
factor e V* 1 ®. Hence the ratio of the probability of being in the first excited state to the 
probability of being in the ground state is 

-3»o/2KT 


" e -flw/2hT 

(b) By the definition of mean value 


= e 


-fiu/KT 


3^-3fiu/2M 


+ 2 

+ e -3h^/2iT~ 


fiu 1 + 3 e ~ ftu !^ gr 

2 " + ’-hu/kr 


6.2 


The mean energy per particle is 

_ He 
e = - 


-dHAr _ 

.-&fcr ' + 7a>r g 


Hence for N particles. 


E = -NpH tanh 


Hi 
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( 1 ) 


6.3 

(a) The probability that a spin is parallel to the field is 

o nH/kT x 


P = 


j-nS/ST + e nH/K 1 “ e -2|iE/kT 


which yields 

For P = .75 and S = 30,CC0 gauss. 


2nH 


k ln[p/l-P] 

T = 3.66°K 

(h) Fran (l), with 4 = 1.4l X 10 -2 ^ ergs/gauss, T = 5*57 X 10 " °K 


6.4 

The power absorbed is proportional to the difference in the number of nuclei in the two levels. 
This is. 


N e' 


jiH/KF 


U e -MH/kT 


n + - n - - e nH/kT + e -v£/XZ ' e nH/kT ~ e “ViH/KT 
where IT is the total number of nuclei. Since pH « kT, we may expand the exponentials and keep 
only the first two terms. 

Hence , 


(i.g -i.jS) 


n + - n_ « N 


HuH 

1 + h5 + 1 _h5 m 


kr 

-1 


kT 


Thus the absorbed power is proportional to T 

6.5 

The volume element of atmosphere shown must be in 
equilibrium under the forces of the pressure and 
gravity. Then if m is the mass per particle, A 
the area, and g the acceleration of gravity, we 
have 


p(z+dz)A 
z + dz 


mn(z)Adzg 


tr 


p(z)A 


p(z+dz) A - p(z)A = dz A = - mn(z) A dzg 


dz 

Substituting the equation of state, p = hKC leads to 
Thus n ( z ) = n(0)e' mez/ ^ r 
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6.6 

(a) 


E 



As T approaches 0 the system tends to the low energy state, while in the limit of high tempera- 
tures all states "became equally probable. The energy goes from the low to the high temperature 
limit when ~ iT. 

(b) The specific heat is, Cy = (^J) , i.e., the slope of E. 



T -* 0, E -> N€ 1# and 


T —* c0 . 


E -> N 



To find the temperature at which E changes from the low to the high temperature limit, we 


evaluate T where 

r 


-(« -€ )/kT 

— H 

E = N €l + | 

e i +e 2 

= N 


2 € 1 

-(€ -6 j/KT 

1 + e 


which yields 


£ 2 _€ 1 


KT 


= In 3 - 1 


or ( e 2~ e l) 88 351 


3k 



The heat capacity is 


. »2 -(« 2 " € 1 )Ap 

r se N <V € i> e 

°V =^T = 


2 “( € 2~ € l)>^^l 

M? Jl+e * x J 


and -» 0 as T -» 0, T -» 00 . 

6.7 

(a) The mean energy per mole is 


E = N. 


where is Avogadro ' s number. 


Oe" 0 ^ + 2^^ 2 N A ee_eAT 


1 + 2e- £ /^ 


1 + to*/** 


( 1 ) 


N. 


(h) We find the molar entropy from the partition function Z. Clearly, Z = £ where 
-e i /kT 

5 = L e , the sum being over states of a single nucleus, 

i 


Thus 


and 


5=1 + e" 6 ^ + e -eAr 


S = k (in Z + PE) = N k In (l + 2e' e ^) + — 


N. ee 


-e/kT 


( 2 ) 


T(l+2e^ 

(c) As T -+ 0, the nucleii go to the ground state, and by the general definition of entropy 

S = N,k In !J = Ul In 1 = 0 
A A 

as T -» °°, all states become equally probable, SI -* 3 and 

S = N A k In 3 

It is easily verified that egression (2) approaches these limits. 

(d) From (l) we have 


c = 


31 a A A ~ 6Al 

oT 


3sT 2 (l+2e" € / Kr ) 
For large T, e € A^ q. and 


c « T 


6.8 



The polarization, P, is defined as N el where 
N is the number of negative ions per unit volume, 
e is the electronic charge, and -t is the vector 
from the negative to the positive ion. Then the 

35 


■e- 


■ X M 



■e- 




o 

l+- 


o 

th 



•eszt x component of P is 


P x = |p| cos 9 = Ne l-t-l -prj- = Ne x 


0f Q O 

vhere x nd 9 are shown in the diagram. The positive ion can have energy + — or , 

and since there are two lattice sites for each energy we have 


eia. 


P 


efa 


P = Ne 
x 


(f)* 2 /+ (~|) e 

e£a e*a 

2 P — P 

e + e 


3 


Nea , . e£a 

— tanh ^ 


6.9 


(a) The electric field is found from Gauss ’s Theorem 

/ E • n da = %q 

J a 

A 

where q is the enclosed charge and n is the vector normal to the surface, a. 
By symmetry, Ea = E 2n r L = %q 

KS) = II £ 

The electric potential is given by 

U(r) = -J E-d l = - |a In 
r o 

Since U(r) = -V at r = R q , we can evaluate q. 

V 


Thus 


E(r) - J In §- £ , 
o 


U(l) 5“ In 


i R 


(b) The energy of an electron at position r is qU(r) = -eU(r) where e is the electronic 
charge. Since the probability of being at r is proportional to the Boltzman factor, we have 
for the density 


p(r) = p(0) exp |^- ^ r A 0 )J = p( °) {jr 


6.10 


(a) The centrifugal force mw 2 r yields the potential energy - 


2 2 
mco r 


-eV/tT ln(R/r ) 


Since the probability that 


a molecule is at r is proportional to the Boltzmann factor, the density is 


p(r) = p(0) exp [mw 2 r 2 /2hT] 


( 1 ) 


(b) Substituting p = N A m > the molecular weight, into (l) and evaluating this expression at r^ 



6.11 

The mean separation, x, is found from the probability that the separation is between x and x+dx, 

_ h£eI 


Where 


P(x)dx 


U(x) = U 


= e 


hT 


dx 


12 61 

© " 2 © 


It is easily verified that the minimum of the potential is at x = a. Since departures from this 
position are small, we may expand U(x) about a to obtain 


U(x) 


- U o + — 


36 U - 252 U - 

- (x-a) = — (x-a) 3 


a a 

where we have kept only the first three terms. It follows that 


P(x)dx = C exp 


36 U (x-a) 2 252 U (x-a) 3 

O ' * O '* 


dx 


a 2 !! a 3 !:! 

We have absorbed the irrelevent constant term, e , into the normalization constant C. 

This constant is evaluated from the usual requirement 


P(x)dx = C / exp 
J 0 


f* 36 U o (x-a) 2 -] p 252 U o (x-a) 3 1 


a 2 kT 


- exp 


a 3 kr 


dx = 1 


The predomi nent factor in the integrand is exp -36 U Q (x-a^/a^krj so the second factor may be 
expanded in a Taylor’s series as in Appendix A. 6. Furthermore, the region of integration can be 


extended to -® since the exponential is negligably small except near a. 


Thus 


1 = C 


exp 


r 36 V x ~ a ) 1 

a 2 !! 


1 + 


252 u (x-a) 3> 


dx 


a 3 kT 


( 1 ) 


where we neglect all but the first two terms. The first integral is evaluated in Appendix A.b 
while the second is 0 since the integrand is an odd function. 

B- 0 . S 

a V 7 r kT/ 

Ey definition, x = f x P(x) dx, and to the same approximation as in (l) we have 


x = 


/ /U \ 5 
— ( 0 \ 
a l 7rkT J 


r 36 u 


exp 


(x-a ) c 


/ 252 U (x-a) 3 \ 

1 + — r n 

K a^kT / 


L a kT 

/ 

The first term in the integrand is a Gaussian times x and the integral is just a. The second 
may be evaluated by making the substitution | = x-a 


3 ? 



X = a + — ; 


U \ 2 252 U 


o 

a\^7rkTy 


3 

a J kr 


h 

S exp 


36 U 


°v 


d | + / a| exp 


G6 u 


O e 2 


ia? 


ae 


Noting that the second integral is 0 since the integrand is an odd function, and evaluating the 
first by Appendix A. we find 

.15 kT N 


x = a 1 + 


U 


Thus 


6.12 


a = 


1 dx 

dT 

x 


.15 ale __ , 15 k 
aU + (.15)akT * U 


since U » kT 
o 


(a) Let the dimensions of the box be x, y, z. Then the volume and area are 

V = xyz A = xy + 2yz + 2xz 

Eliminating x, we have 


( 1 ) 


V = yz 


A - 2yz 
y + 2 z 


( 2 ) 


SV oV 

From the conditions for an extremum, ^ ^ = 0, it follows that 

(A-2yz) (y+2z) - 2yz (y+2z) - y (A-2yz) = 0 
(A-2yz) (y+2z) - 2yz (y+2z) -2z (A-2yz) = 0 

. . A 2 a 2 

Subtracting these equations yields y = 2z. Then substitution in (2) gives y = , z = — r=q , 

A t v3 2 \/ 3 

and from (l) we find x = — . 

V3 

(b) From equation (l), we have 


dV = yz dx + xz dy + xy dz 


(3) 

w 

(5) 


The equation of constraint is 

g(x,y,z) = xy + 2yz + 2xz - A = 0 

and kdg = k(y+ 2 z) dx + \(x+ 2 z) dy + k( 2 x + 2 y) dz = 0 

Adding ( 3 ) and (4) and noting that the differentials are now independent, we have 

yz + \(y+ 2 z) = 0 

xz + k(x+ 2 z) = 0 ( 6 ) 

xy + \( 2 x+ 2 y) = 0 

We multiply the first equation by x, the second by y, and the third by z, and add. 

3 xyz + \( 2 xy + kxz + hyz) = 0 

Substituting (k) yields 

x. = - 2 S£ 

2A 
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and substituting X into equations (6) gives 


1 - g (7 + 2z) = 0 
1 - g (x + 2z) = 0 
1 - § (2x+ 2y) = 0 

A 

The first two eauations show that x = y and the third yields x = Then substitution in the 

a 1/2 a 1/2 k l/2 

first and second equations gives x = , y = , z = — — 

& a/3 2V3 


+ S„ = -k E P ^In P ^ -kZP ^ln 


1 2 


But since E P ^ = E P ^ we may write 


S. + S = -k E E ? ^P ^ 2 ^ln P ^ - kEEP ^ 2 ^P ^ In P ^ 
_l cL r s r sr s 

r s s r 


= -k E E P ^P ^ In P + 
r s r 

T* S l— 


= -kEEP ^P ^In P ^P ^ 


= -k E E P In P 

„ „ rs rs 
r s 


since P ^ p ( 2 ) = p 

r s rs 


From the definition, we have 


S 1 + S, = -k E P ^ln P ^-k E P ^In P 
x d _ r r s s 


Substituting P ^ = E P and P ^ = E P m (x) we have 
r rs s rs v ' 

s r 


S, + S 0 = -k E E ? InP ^-kEEP In ? ^ 
1 2 r s rs r r s rs s 


(1) B (2) 


-k E E P In ? ' ,x ' ? 

rs r s 

r s 


For the total system, the entropy is 


Then from (2) and (3) 


S = -k E E P In P 
rs rs 

r s 


S-(S +S ) = -k E E P in P + k E E P In P ^P ( 2 ) 
'12' rs rs rs r s 


= -k E E p in ( 

r s p ' 'P ' 

r s 
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(fc) Since -In x £ -x + 1 or In x § x-1, we have 

P (Dp (D 


In P ( 2 ) 

S ---' * ^ -w - S - -1 


Fran (!<■), 


rs rs 


P (Dp (2) ( (Dp (2) 

r s skZZp.J 1 1 


_ _ rs \ P 
r s \ rs 


S-(S 1 +S 2 ) = k Z Z P rs In 

r s rs 

But since Z Z P^^P,^ 2 ' = 1, it follows that 

S-( s i* 5 2 ) S k Z Z P^^P^^-l = 0 


/ 


r s 


r s 


r s 


Thus S S S 1 + Sg. 


6.15 


(a) Since S = -k Z P In P and S = -k Z P ^ln P ^ , we have 
\ j r r o r r J 


S-S = k Z 
0 r 


P In P + P In P (D _ P 2n P (D+ P ( 0 )p n P (°) (D 
rrrr rr r r 

where we have added and, subtracted the term k Z P In P ^ . We show that the sum of the last 

r r 
r 

two terms in (l) is zero. 


Since 


P (0) - e"% 
r ' 

In P ^ = -gE - In Z 
r r 


Then using Z P = 1, we have 


-ZP lnp(°)=pZPE + Z ? In 


Z = S3 + In Z 


Similarly 
Adding, we obtain 


Thus (l) becomes 


S-S Q = k Z 
(0) „ (0) r 


(b) Since In 


Z P ^In P -p Z P ^3 - Z P ^ln Z = -?3 - In Z 
r r r r r 

r r r 


-Z? InP^+Z - ^In ? = 

r r r r 

r r 

( 0 ) 


-P In P + P In P 
rrrr 


( 0 ) 


= k Z P In 

r P 
r r 


? r " P r 


1, equation (2) yields 


■S § k Z P 

F ( °) 

— - 1 

= kZ 

p <°>. p‘ 

° r r 

p 

r 

r 

_ r r 


( 2 ) 


= 0 


Consequently S Q S S. 


UO 



CHAPTER 7 


Simple Applications of Statistical Mechanics 


7.1 


,th 


(a) We label the positions and momenta such that r. . and p. . refer to the j molecule of type 
i. There are N. molecules of species i. Then the classical partition function for the mixture 


of ideal gases is 


exp 




(Slrl + ** ,+ 2>-M ) 


3 2 

2m^ 


2 | d3 ^ll'*- d3 ^M k d3 £ll*-- d3 EM k 


^k 


h ^ h 3M * 

h . . . h 

o o 


The integrations over r yield the volume, V, while the p integrals are identical. Since there 


are N^+N^. . .N^ integrations, we have 


3P 


z , = „(v- + Vr r - Sr ai-|Ci + - + V 


htl 


The term in brackets is independent of volume, consequently 

In Z' = (N^+. ..+W k ) In V + In (constant) 


and 


> ~slv^ z ' ' <V— + V k 

pV = (N ] +...+N k )kT = (v 1 +...-*V k ) RT 


(b) For the i td gas, p..V = v^RT, hence by (l) 


p = 2 p. 


(i) 


7-2 


(a) 3y the eauipartition theorem the average value of the kinetic energy of a particle 
€ = | kT. 

(b) The average potential energy is R 

/ mgz e ” mgZ dz 


is 


a -Smgz 


dz 


Then 


5 *'^ to / n 


.-ngz B 3 B - 1 

8 ** • "Se ^ — = 5 ^ — 


On carrying out the differentiation, we find 


u = kT + . — 

2_ _ e ^gL/iT 


1-1 



7.3 


(a) Before the partition is removed, we have on the left, pV = vRT. After removal the pressure 
1S -n - RT _ 2g_ 

P f (l+h)v 1+b 

(h) The initial and final entropies of the system for different gases are 

S, = vR |!n ^ + | in T+a 1| 


+ vR 

In — — + ^ In T+o- 


N. v 2 2 

A _ 


vR | 

In + | In T+a 

+ vR 

in ASH + 3 J 


N. v 2 1 

A 


N a v 2 2 


Here one adds the entropies of the gases in the left and right compartments for while S f is 
the entropy of two different gases in volume (l+b)v. 

Then 


AS = S -S. = vR 
f i 


Kiiil . u JL . ^ w. 

V V V 


vR In 




(c) In the case of identical gases, is again the sum of the entropies of the left and right 
compartments. is the entropy of 2v moles in a volume (l+b)v. 


= vR 


S f = 2vR 


V R 

In - — + ^ In T+a 
N^v 2 o 


vR 


bV R 

In + — In T + a 
N A V 2 ° 


m + j T + 

2N^v 2 


Thus 


AS = vR 


In 


( 1+b ) v .1.X.1.W 

2» A » B a v B a v 


vR In 


Sb 


7.5 

(a) The system is isolated so its total energy is constant, and since the energy of an ideal gas 


depends only on temperature, we have 

AS X + A3 2 = CyCT^) + C y (T f -T 2 ) = 0 


or 


T = 

f 


T 1 + T 2 


The total volume is found from the equation of state 

V i m i V 2 HE 2 

V = — — — + — — - 


Thus the final pressure is p^ = 


( V]L +v 2 ) RT f (v x +v 2 ) 

V " = 2 U v 1 Ti / Pi)+( v 2 T 2 7p 2 ) 


T l +T 2 
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y 


(b) Using ^ •= ~ , ve have for the initial and final entropies of different gases 



kT , 



s i = V 1 R 

In + 7 In T. + a. 

P-L 2 11 

+ VgR 

In + In T„ + c_ 

P 2 X 2 2 


S f = v x R 


k / v l^l v 2^2\ 3 ^ 1^2 

V A p i % ' 2 2 11 


+ v 2 R 


„ k / V 1 T 1 v 2 T 2\ 3 , T l" Hr 2 

““d — + — ) + a^ — + °a 


Thus AS = S f -S i = v x R 


V 2 T 2 P l\ 3 T 1 4T 2 

111 ,1+ -h 


+v 2 R 


v 1 T iP 2 \ 3 T ±^o 

In 1 1+ ■ | In -i- 2 

v 2 T 2 P 1 ' 2 211 


(c) For identical gases 


= v,R 

kTi ^ 

In + — In T.+a 

+ v_R 

kT .. ^ 

In + ^ In T„ + a 

1 

P x 2 1 o 

2 

P 2 2 2 o 


= (Vj+v^R 


In 


V 1 T 1 V 2 T 2 \ 3 T 1^2 

+ ■—■) + % In + a 

^ o 


AS = v-jR 


(Vi^v 2 ) l 1 T-jP, 


, 2/ 2 2^ 


(VV V p l p 2 / 2 

+ VgR 


In 


v 1 T 1 P 2 \ 3 T i ^ 2 
) 4 ln ± 


W V 2 Vl / 2 2T 2 


7-5 

We take the zero of potential energy so that if a segment is oriented parallel to the vertical 
it contributes energy Wa and if antiparallel it contributes -Wa to the total energy of the rubber 
band (Thus if the rubber band were fully extended, the total energy would be -RWa). Since the 
segments are non-interacting, 


-Wap .^WaS 

= N - 77 . ' - = Wa tanh 


-Wap + e WaS 


kr 


7.6 


Since the total energy is additive, 

F 2 

V Mp.) + u (V C -n) = + u (V°-n) 

where U is the energy of interaction, the eauipartition theorem still applies and 


€ = | 


7-7 


If the gas is ideal, its mean energy per particle is 


2 2 
T3 p 


and the mean energy per mole becomes 2 = W^kT 



7*8 

For clas sical harmonic motion the mean energy is 

*-£♦ t- 2 ? 

5/ equipartition, each of these terms contributes ^ hr. 
Hence p °= x « T 


7-9 

(a) 

The mean elongation is found by equating the gravitational and restor ing forces 

ax = Mg, thus x = 

a 


(b) By the equipartition theorem 

1 / — v 2 

p a (x-x) 


| KT 


or 


(x-x)' 


kT 

a 


( c ) 


7-10 


(x-x) 2 


2 _ 

= X 


m 2 or M = 


.a, 


a 


(a) Let the restoring force be -ax. Then the mean energy of N particles is 

— 1 »2 1 p 

E = N (r mx +^ax ) 


By equipartition 


E=If(|kT+|kl)=NkT 


c =g = Hk 


(b) If the restoring force is -ox3, the mean energy per particle 


is 


- 1-21 T 

e = mx +^ax 


The kinetic energy term contributes — kT by equipartition. The mean of the potential energy is 


found from 


ax 

T 


4 -, 


exp 


exp 


Bax 

-rr 


ax 

T 


dx 


Bax 


47 


dx 


44 



T 

ax 

T* 




la 


r e -p(«» te = .ur 

-03 ^ ^ “03 


e -<^ A p'lA fly 


where we have made the substitution* y = p 7 x. 


iA 


Thus 


T 

coc 

T 


¥ 


(- i In p + In 


s -oy A 


dy} 


= ¥ K 


Hence 

Then the total energy is E = He and 
7-11 


_ 1 -2 1 — E 1 1 S 

e = imx +Jcoc = |kT+£kT=|kT 


C =|| = t m 


Since the probability of excitation is proportional to e 


nfito 

IT 


where n is an integer, it is 


clear that the two parallel components are negligibly excited at 300 where » 300 k. 

The perpendicular component is excited since « 300 kj and since 300° is in the classical 
region, the equipartition theorem holds and we have for the mean energy per mole 

E = N a (| mi 2 + | mUj 2 x 2 ) = H A kT 


Thus 


7.12 


0 * § * V ■ 5 


(a) Consider a cube of side a. The force necessary to decrease the length of a side by Aa is 

Q O 

k Aa and therefore the pressure is Ap = k Aa/a . The change in volume is AV = -a^Aa. 
o o 


Hence 


1 f AV 

* ~ " V Up. 


2 . 

a Aa 


k Aa/a 


a_ 

K 


( 1 ) 


(b) The Einstein temperature is 0 E = ftco/ k. Since to = where m is mass, we have from (l) 

( 2 ) 


x K „ 

s «*.._£ = t f_ 

E k i m k Vmn 


The length a is found by considering that the volume per atom is n/pU A where pi is the atomic 


weight, p the density, and N A is Avogadro's number 


Thus 


1/3 

a = (~ J ~) = ( 


(8.9)(6xl0 23 ) 


1/3 Q 

) = 2.3 X 10 cm 


Substitution of m = pi/N and the given values yields 


0 E » 150 K 


^5 



for J = — "becomes 


7-13 

Ve have 2^ = N Q g (i Q J Bj(tj). By equation 7.8.14, ^ (tj) 


Si(t) = 2 (coth tj - j coth |) 


cosh n _ cosh (n/2) 
sinh 7] sinh (t]/2j 


= 2 


Thus 


7.14 


2 2 

cosh (n/2) + sinh Ct/2) cosh (b/2) 

2 sinh(r)/2) cosh (q/2) " sinh (t)/2J 

IliT^fcosh (n/2) ' = c ' o ' f {^2] = tanh (n/2) 

N 0 S u o . , s ^o H 

\ = 2 


The energy of the magnetic moment in the field E is 

2 = -ji*H = -|iH cos 9 

where 9 is the angle "between the magnetic moment and the direction of the field or 2 -axis, 
the probability that the magnetic moment lies in the range 9 to 0 + d0 is proportional to 
Boltzmann factor and the solid angle 2 jt sin 9 dQ. 


Thus 


P(e)d 0 « sin q 


and 


'7T 


IJ 


cos 9 sin g cos 0) 


«Z ■ Vz * 


r 

n 


BliH cos 9 . . 

e Kh ^ sxn 9 dS 


or 


N 


,, o S . f B(iH COS0 

“z = T ¥ W e 


N .. BpH -3uH 

sinf5 de = f Ip 111 “ — FF 


7-15 
We have 


0 5_ ^ 2 sinh BuH o . 

E eg puE _ hb sinh p|iH 


pEgcosh 3pE -sinh guE 


“z = V 


coth PuH - 




«z ■ R o 8 1*0 J 


U) 


46 


Then 

the 



lere By (7-8. I*)-) 


B J^ J 


(J + §) coth (J + hn - | coth | r, 


If T] « 1 and J » 1 in such a way that Jt] » 1, B (f]) Becomes 


Bj(t)) = J J coth jT l " \ (|) = coth J, 1 ' 


where we have used 


coth x = 


x , -x 
e + e 

x -x 
e - e 


+ -.. ' • - ! . •? + (1 - x + w 1 for small 3 

(1 + X + ...j - Xl “ X + ...) X 


Let pTTp = Jtj = Jga Q Hp where jj. = gh Q J is By (7*8.2) the classical magnetic moment. 


Becomes 


7.16 


*Z = V 


coth ppH 


Pul 


(a) The energy of a magnetic moment in a field is E = -4 • H. For spin — atoms 

fc <V H i> 




uH 2 P 


dH 1 P 


= exp 


(*) 

n(Z 2 )/n(Z 1 ) > 1 since Hg > ^ , 
(c) If uH « hT ve have 


n(Z_) ^ 

v 2.' _ e +_e 

n^) ~ pH-B -fjH 7 P 

e + e ~ 


~ 1 + fc (h 2 _h i ) 


and 


r-(Z 2 ) 

^(zj 


(l + pHgB + i lA 2 V) + (l - hEgP + ^ U^p 2 ? ) 
(1 + fiHjp + | p 2 ^ 2 ? 2 ) + (1 - nE^ + | u^V) 


2 2 ? 
2 + u Hg P 

“ 2 _ 2.2 

2 + n p 


1 + l <et> 


1 

1 - - ( — -) 

2 v hT ' 


n(Z 2 ) 


( z^T “ 1 + 2 " H ! ) 


7.17 


To find the fraction, of molecules with x component of velocity Between -v and v. 
integrate the distribution Between these limits, i.e.. 


47 


Then (l) 


we must 



i r v 

| = = / g(v )dv = 

n J ~ ® v x' x 


m 

/2kT 

m 


^27rkT ; 


■(“v /2KT) 
e dv 


Making the change of variable, y = ,^£L v , we have. 




5 = f 3" (y dy = f e"( y / 2 ) = 2 erf VF 


VFr J-y/2 -/irr ^ 0 

Then from tables | = .8^3. 

7.18 

In problem 5 *9 w e found that the velocity of sound is 

- = 


where 7 = C and 4 is the atomic weight. Since 4 = N^m,we have 


u = (^) 2 = ( 2 * S ) 

v N.m ; ' m ‘ 

A 

The most probable speed is v = (2kT/m) 2 . Thus 


u = (f ) 4 v 


For helium, 7 = 1.66 so that u = .91 v } and the fraction of molecules with speeds less than u is 




.9iv 


3 2 

mv 


0 


F(v)dv = % 


.9l(2M/») 2 2 2 - m 

<SB> T 8 


dv 


Making the change of variable y = (ra/kT) 2 v, we have 

% 


5 = 


(27T) 


^ 2 -y 2 /2 fl 

y e ^ ' dy 


'0 


f a -y 2 /2 

This integral may be evaluated by noticing that integration by parts of / e J ' dy yields 


0 


£ e' y2 / 2 dy = ye“ y2 / 2 


+ / 

0 J 0 


2 -y /2 . 
y e J / dy 


Since the integral on the right is the required one, we have. 


.1 -J r .91^ 2/ 1 - ■$ r .91V2 2/ _ P -1 

I = % 2 (2) 2 J y 2 e " y / 2 dy = % 2 (2) 2 J e 7 / 2 dy - .91V2 exp|-(.91 2)2/2) 


1 - I r * 91^2 


The integral on the right may be evaluated from tables of error functions. Thus we find | = .37* 


h8 
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(a) 


v x = 0 "by symmetry 

(b) v 2 = — by eauipartition 

v ' x m 

(e) (v 2 v } = (v ^ + v 2 v + v 2 v ) = 0 
v ' '•X'' x y x z x' 


(d) (v^v)=v' 5 v =0 

' ' X y' x y 

(e) (v +bv) 2 =v 2 +2bv v + b 2 v 2 

w x y' x x y y 

= (i + b 2 ) 3£ 

2 m 

( f ) v 2 v 2 = (^) 
v ' x y v m 1 


o /o mv 

(a) ^ ^ dv = J f % (^) ve ^ dv 

l. l. 

By appendix A. 4 ve find (1/v) = (^-) , and since v = (— — ) , 

VliCl 71 EL 

^ = 1.2J 

dA) w 

12 2 

(b) Since the energy is e = ^ w , ve have v = 2e/m and dv = de/ V2me. On substituting these 
relations into the speed distribution F(v) dv, it follows that 

_ 3 i _ 

F(e)de = 2?m (7 tKT) 2 e 2 e M de 


The most probable energy is given by the condition = 0 or from problem 7.20b 


e_ jl 
1 4 a ICT _ e 2 


e M =0 


The most probable speed is 


Hence 


~ 1 

e = | 1ST 

~ 1 

v = (2kT/m) 2 

1 ~2 

- nrv = kT 


V = V (l + — ) = V 
o c o 


since v =0 
x 


(b) The rms frequency shift is given by 


(Av) ms = _ (v_v ) 2 _ 2 = / ‘ y2 _ 


—x _ / 2v v 

v - v o + ^ 


Since v^ = 0 and v^ = kT/m by ecuipartition, we have 


2 2 /. kT 

v =v o 1+ ~ 
\ me 



/ . \ r 2 W 2 il 2 

(Av) = + — 5 - V - V 

rms o 2 o o 


r o /iT' 


(c) The intensity distribution is proportional to the distribution of the x component of 


velocity. Since v = — (v-v ), we have 
X v 'o ' 3 
o 


2 , >2 

r me (v-v q J - 

l(v)dv = I exp — dv 

L 2kT v J 
o 


where I is the intensity at v=v . 
o o 


(a) The number of molecules which leave the source slit per second is 


<J> A = r 1 nvA = 
o 4 


= 1.1 X 10 molecules/sec 


where A is the area of the silt. 

(b) Approximating the slit as a point source, we have by ( 7 . 11. 7) > the number of molecules 
with speed in the range between v and v + dv which emerge into solid angle d2 is 

A$ (v) d^v = A[f(v)v cos 0] [v 2 dv dfi] 

cos 0 = 1 for molecules arriving at the detector slit; hence the total number which reach the 


detector is 


« r- & ^ 


A fl d p s 
7T V 2nmKF 


where we have used the results of appendix A,k and p g = nkT. Letting L = Im be the distance 
between source and detector, the solid angle is approximately A/L . Substitution of the 
given values into (l) then gives 1.7 X 10 11 molecules/sec arriving at the detector slit. 

(c) In the steady state, the number of molecules entering the detection chamber per second is 

equal to the number being emitted. Prom (l) we find 

,2 . - 

A Pq A P., 

~ S ■ ■ = 0 d A =-— 4r- 

7rL y27nnlfi? i/Stt mkT 


_ A 

P d _ P S ,.2 


2.4 X 10 inm of Hg. 


The rate of change of the number of particles in the container is given by 


f = - cpA = - J nvA 


50 



Since the pressure is proportional, to the number of particles, we have 


dp vA 
dt = " W P 


p(t) = p( 0 ) exp 


vA 

tw!| 


where p(0) is the original pressure. Hence the time required for the pressure to decrease to 


p(t)/p( 0 ) = l/e is 


t -12 


As in problem 7*2^, we have 


dn vA 

dt = " 5v p 


4- W , P 

t In — 


where p(0) is the original pressure. The mean speed, v, is for nitrogen at 300°K, 

v = JPf = .^75 x 10 5 cm/sec 


.... to Vil . If . 07 sec 
(.^75 x lCT)(l) 10 _i 


(a) Since the rate of effusion is <p = p/(27r m kT) 2 , and the concentration is proportional to the 


pressure, it follows that 


v = f2 (5T 

c l* ^ V 


(b) For the uranium separation, we have 


A. 

c '233 // ' C '23Q = C 235 //C 238 ^235 + o[l 9 ^ = ( c o^/ C o-3ft) (1*003) 


23 V 238 ; 


The rate of change of the number of particles inside the container is 


dU 1 - NA / 8 kT AN 


dt t 


nvA = 


77 m Vm 


thus defining A. Since pressure is proportional to the number of particles, we find after 


integrating 


/ At 

p/p o - exp b ®J 


For Helium gas p/p Q = l/2 at t = 1 hour. Substituting we have 


A 1^2 
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where m^ is the mass of the helium molecule. From (l) and (2) it follows that the ratio of the 
Ne and He concentrations in terns of the initial concentrations is 

■vV - < W 0 [- x • 1/ vv') ^ 

* <v/V> 0 exp [-^ 2 (V"He/v’ *] 

Since (n^/n^) = 1, we have at t = 1 hour. 


^Ke/^He ~ 2 


( 1_ Ie^ ^ ( 1_ A AwSfe ^ 


= 2 


7.28 


(a) The rate of change of the number of particles in the left side of the box is 

® L (t) 


dt 


= (no. entering/sec) -(no. escaping/sec) 


. 1 


Since the rate of effusion is ^ nvA, we have 

dN x (t) 


dt 


= W/2 [ N 2 ^)-V t ] = § ^^(t) 


where N is the total number of particles. From the equation of state, 

Ni(-t) = P x (t) ^ and N = 


V „ P l^ 0 ^ + P 2^ 0 ^ V 


kT 


(1) 


Thus 


dp, (t) 

dt = 2V ([fl^ + P 2^ “ 


^(0) + p 2 (0j^ A(0) - p 2 (0)\ - f t 


Integrating,we find p^(t) = 

(b) The initial entropy is 


S. . ^(0)* (in + | In T + a^) + B 2 (0)k 


in 


^15) + f ^ T + °o) 


In final equilibrium, = N 2 = 2 

S f = m (In 5 + | In T + c o ) 

Then AS = S^-S^ and using (l), we find 


Pn (0) + Po(®) y 
AS = -± — -£ 1 In 


2kJ 

p 1 (o)+p 2 (o) 


p l(°) v 

sr In 


2 T 


KT 

P>) 


P 2^°) Y - 
2 T 111 


~ kT 1 


V f „ 3?i(o) 2P 

21 (p^o) m Pi ( 0) + p 2 ( 0 ) + p 2 ( °) 111 i^ToJ 


2Po(0) 


+ Pp(o)y 
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(a) Inside the container = 0 by symmetry. 

(h) The velocity distribution, 0(v), of the molecules which have effused into the vacuum 

3>(v) d\ = f(v) v^ d\ 

where f(v) is the Maxell distribution. 2 


v = 
z 


' dv / 


' dV xJ 

■<0 -<o 


dv / dv $ (v) v 


dv / dv $ (v) 


, 2kT 2 

dv e v 

z z 


PVT 

dv e v 

z z 


The second equality follows because the integrations over v^ and v^_ are identical in the numera- 
tor and deno m i n ator and therefore cancel. The integrals over v are tabulated in apnendix A A. 

z 


1 /2kT\ 

2 U/ 


ZE 

(2 m 


(a) In time dt, molecules with velocity component v which are in a cylinder of height v dt 

z z 

with base at the aperture can escape from the container. Molecules of higher speed can escape 
from a larger cylinder than can molecules of lower speed, and a proportionally larger share of 
the fast molecules effuse. Thus the mean kinetic energy of particles in the beam is greater 


than in the container. 


1 2 

£q = g mv = 


'3 1 2 / \ 

d v g mv $ (v; 
d^v 0 (v) 


o 2 

Since <&(v) = f(v) v cos 9 and d J v = v dv sin 9 d0 dcp, we have 

r® rV 2 n't To . i 


P 00 P 7T/£ P7T ^ i 2 ?Vp 

/ / / (v dv sin 6 d9 dcp) mv ) e v cos 9 

j o ^ O J o 

— my 2 

r°° r 71 / 2 of 2 " pm 

/ / / (v^ dv sin e d0 dcp) e v cos 9 


1 f 5 2kT . 

r i v e dv 


= 2kT 


TV 

^ o 


2kT , 
e dv 


where the integrals over v are evaluated by appendix A A. Since the mean kinetic energy in the 


53 



the enclosure is e. = 7 : kT, it follows that e = ^ e. . 

i 2 o 3 i 

7.31 

When a molecule is scattered, it commutes momentum Ap = 2mv to the disk. The force is 

z z 

Ap / dfc. To find the total force, we must multiply Ap /dt hy the number of molecules with 
z z 

velocity in the range (v, dv} which strike the disk in time dt and then integrate over all 
velocities. Clearly only molecules in the cone of half -angle 0 q = sin ~ r/l will strike the 


disk. 


F = / d J v 


A$ (v) dt 


1 A P Z 
dt 


-X 


“ rS o r 2 *, 2 


(v av sin 0 d0 dtp) (A f(v) v cos 0 dt)(- 


o u o 


2m v cos 0 X 
dt > 


(2mA) (2rr) (- j cos^0) 


0 . 


v f (v)dv 


o u o 

The integral can he expressed in terms of the mean square speed inside the container, i.e.. 


v^f (v)dv = j£- 


mx. „ 1 2 . , 1 3„\ 0 1 2 A , 3,-1 Rx 

Thus F = =■ nm v A (- w cos 0) = — nm v A 1 - cos (sin ■=■) 

3 3 o3 ' L 


But hy eauipartition we have 
Hence 


1 2 . _ — 
— n mv = nkT = p 


F = pA 


3 , . -1 R\ 
1 - cos [sm — ) 


CHAPTER 8 

Equilibrium between Phases and Chemical Species 


8.1 


By equation (8.3.8), we have the proportionality 


CP(v,T) dV dT « exp -G (V,T)/kT 


dV dT 


where V and T are the volume and temperature of the portion of substance of mass M, and T q is 
the temperature of the surrounding medium. Expanding the Gibbs Free Energy G q (V,t) about the 


minimum at V = V, and T = T yields 

Bg 


Bg 


b 2 g_ 


G q (V,T) = G n (V,T) + (sr) (T-T) + («2) (V-V) + | ( 1) (T-T) 2 


v bt 

o 2 G 


BV 


BT 


o 2 G_ 


( 1 ) 


+ I ( ^2 £) + (aTS?) (v-v) (t-t) + ... 


BV BT 7 


5 ^ 



where the derivatives are evaluated at T = I and V = V. Since 
the first derivatives in (l) are 0. By (8.1;. 7 ) 


ve are expanding about a minimum. 


Since this is 0 at T = T, we have T = T q . (Se/oT)^. is the heat capacity, C . 


o G \ T 


T SC.. 
0\ , V\ 


3 ' 3 <V * & - fH-g 


and at T = T = T, 
o ' 


From ( 2 ) we also find the mixed derivative 


2 T 


T SC T , 


(fv>„ <sr>„ - n - rKsrl - 0 


T V 


at T = T. 


By the equation preceding (8.1;.13) 

cG 


(^) = (T-T )(^) - P + p 

cV ' o y 'cV-' f p o 


- (M « ) 0 


-(f) = “ (f ) at T = T , V = V 

u T T 


1 SV 1 

By definition, k = - — (2=) , thus ( 5 - ) = - — . Substitution of these results into (l) yields 

v Sp T \svy V k 


G o (V,T) = G (V,T ) + (T-T 0 ) 2 + (V-V ) 2 

~o 2 Vk 


and since V = M/p Q and = Mc^ where c^. is the specific heat per gram, it follows that 


(?(V,T)dV dT °= exp - 


rs (IJr / ' asrw < T ' w /°o > 2 ® ® 


where we ha ve absorbed the constant term G (V, T ) into the normalization factor. Performing the 

o o 

usual normalization, we find 

/cl p \ 2 Me p p p 

<?(v,T) dv dT = exp "^ (T_T o) ■ ^o* dvdT 


(a) For the solid. In p = 23.03 - 375 V T and for the l ia - uid In P = 19* ^9 - 3 C 63 /T. At the 
triple point, the pressures of the solid and liquid are equal. 

23.03 - 375 =19.^9 - 306 3/T 
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hence 


T = 691/3. 5h = 195 °K 


y 

(h) Since In p = - — + constant where •£, is the latent heat, we have 
RT 

^sublimation = 3 ^ R = 31 ’ 220 j^ es / nole 
^vaporization 3 3063 R = 25 ^ 8 ° j oules / mole 

"^fusion ~ ^sublimation ^vaporization 

= 31,220 - 25 } b&Q = 5740 joules/mole. 


The latent heat of sub lima tion is equal to the sum of the latent heats of melting and vapori- 
zation, i.e., = ■£, + Since -t = TAV (dp/dT), the Clausius -Clapeyron equation, where v is 

the volume per mole, it follows that 

f— = 7 . T (v - v ) + T fv - v ) f—®’) 

'dT' TTv - v J o '‘lie sol A dT ; o v gas lia A dT' 
s o v gas sol' L m * vj 


= f 2. _ 

'dT' Ip p 

s \ o K s. 


T fv - v ) + T fv - v ) 

1 o k lia sol' 'dT' o v gas lia n dT' 

_ * m v 


Since ^ ^ and AS -» 0 as T -» 0 "by the third law, it follows that 0 as T -* 0. 

dT AV dT 


(a) Q/L moles/sec evaporate and are swept out "by the pump. Then "by the equation of state 


£ V J±L 

L RT 


or V TT 


(h) We find the temperature at the pressure p m from the Clausius -Clapeyron equation 

dr _ L _ L 

dT _ t(v - v. . . . , ) ~ T v 

' gas liquid' gas 

where we have neglected liquid the volume per mole of the gas is much greater than that 

of the liquid. Then using v_ &s = RT/p we find 

r 5 m dp L / T m dT 


1 T T - 

o 


*m L /I 1 > 

^ p R T " T 

o o m 


Hence 


T = T ( 
m o' 
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8.6 


The intensity, I, is given by 


1 — p 

I = £ nv = 


V2 tt m kT 

Differentiating -with respect to T and dividing by I, we have 

1 dl = 1 d£ 1_ 

I dT p dT “ 2T 

Since p is the vapor pressure, the Clausius -Clapeyron equation yields dp/dT. 

L 


(1) 


dp _ ^ L 

dT ~ T(v 


v~. ZJ T v 
gas liquid' gas 


Substituting into (l) and using V g as = HT/p gives 

1 dl _ L 1_ 

I dT ^2 “ 2T 

8.7 

Differentiating the molar latent heat yields 

41 - <#) § - (i> 


dT 


cip' dT 'ST' 

3? 


(i) 


Now -L = t(s 2 ~s^) where (s^-s^) is the molar entropy difference between the phases. 


Thus 


<i> ’ 


8s 8s 

^5p"^ " 

T J 


But (os/8p) T = -(8 v/8t) , a Maxwell relation, and since the coefficient of expansion is 
a = -l/v (cv/ST)^ ,it follows that 


( ll ‘ ^ [“2 V 2 ' Vl 


( 2 ) 


Differentiating at constant p gives 

# - <V S 1> + 1 


CS_ CS- 

■ ( §r ) 


T(8s/oT) is the specific heat at constant pressure. Thus 


<#> -&+(< 


) 


■2 *1 


( 3 ) 


Substituting (2) and (3) into (l) yields 

= -T(a 2 v 2 - o 1 v 1 ) § + % + ( c p 2 - c p i ) 


dT 


Since 


dp _ -t 
dT ~ T(v 2 -Vi) 


dt , x . I (a 2 V 2 ‘ ¥l^ 

dT v P 2 p x ' T v 2 -v 1 
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3.8 


The bar moves "because the higher pressure underneath it lowers the melting point. The water is 
then forced up along the sides to the top of the bar and refreezes under the lower pressure. 

To melt a mass dm of ice, heat -Idn must be conducted through the bar. The heat passing through 
it per unit time is AfT1 

51 = - K (be) f 

where AT is the temperature drop across the bar and k is the thermal conductivity. Letting x be 
the position in the vertical direction, we have dm = p^(bc dx), and equating the rate at which 
heat is conducted to the rate it is absorbed in the melting process we have 

-K (be) f=l P . (bC g ) 

AT is given by the Claus ius-Clapeyron equation 


(1) 


dr __ Ap _ 
at ~ AT t(v 


l 


- v. 7 


water ice' 

where v is the volume per gram or l/p. Then since Ap = 2 mg/bc. 


Substitution in (l) yields 


8.9 


-6 be 'p p. 


dx _ 2 mg kT / I 1 \ 

dt ~ . p2 ( ‘p. " p ‘ 

aoc 'Op. i w 


From the Clausius -Clapeyron relation, dp/dT = L/TAv we find 


r dT _ f AV 

J T'J L 


dp 


We know L(cp) and AV(cp), and we can find dp/dcp from the curve of vapor pressure vs. emf. 


Thus 


273.16 


dl* 

T' 


AVfaO ,dpv 
L^cp*) Me?' 


dep* 


where 273 «l 6 °K is the temperature of the reference junction of the thermocouple. 


Then 

8.10 


T = 273.16 exp 


■AV 

L 






We divide the substance into small volume elements. Tie total entropy is just the sum of the 
entropies of each of these, or 

S=ZS. (E ± ,N.) 
i 

where E_. and denote the energy and number of particles in the i^ volume. If we consider 
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.th 


an interchange of energy "between the i 0,u and volumes, we have since the entropy is stations 


3 S. cS - 

dS = 0 = vtF dS. + dS . 

0 


diii\ 2. 


"but since dE. = -dE.« ve have 
i a’ 


cS. oS. 

oEi ^ 


or 


1 . = T. 


Similarly, if we interchange particles 


implies 


8.11 


Ss. as. 

1 = __2 
SnT oiT 

i j 


cs. as. 

dS = 0 = v^r 1 dN. + ^ dli. 

oil. a oH. j 


or 


^i =ti j 


(a) The chemical potential can "be found from the free energy by n = ^ . We consider a small 

oN 

volume element at height z over which the potential energy may he considered constant. It is 

clear that the addition of a constant term to the energy does not change the entropy of the 

substance in the small volume element. Explicitly, the partition function will he the free 

— 3 7 np r z 

particle function times the factor e . 

-pmgz^ 


Thus 


S = 


k In Z„ e 
free 


+ 3 l E free + ***' 


= k 


In Z_ + 3E„ 
free free 


Hence the free energy is, letting N he the number of particles in Ay 

F = S - TS 


ana 


ip 5 

p = ^ = ■! hr + mgz - kT 


= ^ llkT + Ifingz - HkT 


vp -s 

In — + ^ In T + a 


iln f + | m T + c c 


where we have put AV/u = kr/p. 

(h) Since p is independent of z. 


dp KT an 

— = mg + — t=- 
az p az 


or 


2£= - ^dz 
n kT 


Integration yields 


mgz 

IT 


p(z) = p(C) e 

8.12 

The relative concentration remains the same. The dissociation is described by 

n C0 ^H-O 
2- = X (T) 

n cc 2 X n 

Generally, when the number of particles is conserved, equilibrium is independent of vc 
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8.13 


The intensity of the I atoms in the "beam is proportional to their concentration in. in the oven. 
Since dissociation proceeds according to 


We have 
which gives 

by the equation of state, 
into (l) 


I 2 = I + I 

a l'\ ’ 

*l\ ' TS^ * V T > « 

The total pressure is p = p + p which yields on substitution 

X 2 

P - Pi = Kp Pi 2 


or 


P z = (-1 +*Jl + 4-iy?' )/2K p 


Thus, if the total pressure is doubled, the ratio of final to initial intensities is 


Of (p I ) f -1 +yi+8K~p 

°i -1 +Vl + lt-K p p 

Since K is small we may expand the square roots to obtain 
P 

4> f -1 + 1 + 4KpP 

~ -1 + 1 + 2K n = 2 

i p* 


For chemical equilibrium we have 


-1 “2 — " K J T ) 


where n is the number of molecules per unit volume. Since n^ = p^./kT -"k follows that 


b. b b 

12 m 

!?/■« i • • 


"b_ +b + • • • +"b 

= (SE) 1 * a K n (T) = KjT) 


3.15 

3y the first law aO, = dE + p dV. Letting d(pv) = V dp + p dY we have 

■d Q = dE + d(pv)-v dp 
= d(S + pV) - Vdp 
= dH - V dp 

where the last follows from the definition of enthalpy, H = S + pV. Then if p is constant 

d Q = dH 


.e., in a constant pressure process, the heat absorbed is equal to the change in enthalpy. In 


a chemical reaction, Q = LE = Z bJn since enthalpy is an extensive quantity. 
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8.l6 


~3y the equation preceding (8.10.25) 

Zb. In £. * = Zb. in n. 
i 1 1 i 1 

Substitution of n^ = p.j/fcT into (l) yields 

Zb. In p. = Z (b. In t. » + b. In hr) 

. i i . ; i s x x ' 

1 n ' 


( 1 ) 


or 


■where 


D D D 

P 1 p 2 •"P m ”=V T ) 

In K^(T) = Z (b d In £ « + b In kT) 


b. 


Thus 


ffl to K.W - s (\ 5 to Jj' + -X) 


e. h. 

By (8.10.30) -=• In £ . * = where e. is the mean energy per molecule. Then multiplying ~ 

CLL 1 J 5 .J. 1 T 


p.V. 

^ 1 = TW vs have 

i 


p.V. 

-1 * 1 . 


b.h. » 


h **,<*) if). 


where In 1 is the enthalpy per molecule. Multiplying and dividing by Avogadro’s number gives 

h * V T > = E ^ ^ 


2 ~ 2 
i ET ET 


where lu is the enthalpy per mole. 

8.17 

We let | be the fraction of R^O molecules which are dissociated. Since the pressure is pro- 


portional to the number of molecules, we have 




2p„ 


s = 


P H p 0 +p H p P H„0 +P Hp 


( 1 ) 


The total pressure is the sum of the partial pressures. 


P = Ptt + Pn + P- 


H 2 t ^0 2 • %0 


Substituting (l) and ( 2 ) into 


2 

p O 

^(T) = -2—Z 
p a,o 


( 2 ) 


we find 


V T) = TTIv 


.3 


(i-s) (2+6) 


( 3 ) 
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8.18 


In problem 8 .l 6 it was shown that 


AH = RT 2 ^ In K p (T) 


Using Kp(T) from problem 8.17 would, give AH for the dissociation of 2 moles of HgO since the 
reaction is 2HpO 2Hg + 0^. Hence for the dissociation of 1 mole 

** = h ^ 

where K^(t) is given by equation ( 3 ) in 8 . 17 . Since f « 1, we have 


V T > - 


Pi' 


(1-0 (2+0 


pg ~ 
2 


Hence 


m ~ ~ dl 111 ~ " 2 dl 111 5 


^ In | is found by plotting In £ vs. temperature using the given values and measuring the slope 
of this curve at 1700°. We find In 5 = .CC 6 h. Thus at 1700°K 


AH = | (1.99)(1700) 2 (.OC60 = 56,000 cal/mole 


8.19 


►N 


(a) The partition function is Z = where £ is the partition function for a single molecule. 

l\l • 

For a liquid we assume a constant interaction potential -tj and that the molecule is free to move 


in a volume N 


v 


Thus 


!, = 


’a 3 „ a 3 r 


where the integration is over Ngv Q and all p. 


s = Vo e “3 f d3 P e ’ P ~ /2m = \ ^ V 


thus defining £ ’ , 
Hence 


No 




= nTT = nTT h e v o ^ 


(b) For an ideal gas, the partition function is 


% 


N 

z s ■ d M 6 


(1) 


The Helmholz free energy is 

F = -KT hi Z = -hi (in 


N 

[V] S 


- In IT J 


= -3sT (N In V t’-N InN +N 

Vs g b g g 
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where ve have used Sterling’s approximation In J ! » u InN - w 

g g g g 


Then 


SF 

= Sir = -lsr v g 5* “ ^ N g ) 

& 


(c) For the liquid, ve have by (l) 

*£,= - 1 


(In 1 

To e &T) v t’ 

V L 

'V o _ 


IT 


X 


- In No 




= -IdT (H£ In v^' + pT) + Ng) 


sf 


£ 


^ = ^ = -kT (ln v 0 S' + (*l + 1) 


Therefore 

(d) When the liquid is in equilibrium with the vapor, the chemical potentials are equal or 

-10? (in V f ’ - In N ) = -KT (in v £’ + prj + l) 
o & O 


Thus 


but since p/kT = N /V ve have 
g g 


li. - v oP 7 ! 4-1 

N o “ 

6 


P _ e 
IT ev 




( 2 ) 


(e) We calculate the molar entropy difference from the Clausius -Clapeyron relation. 

dp As As 

dT ~ Av ~ v 

g 

where in the second equality ve have assumed that the molar volume of the liquid is much less 

than that of the gas. Since pv = ST ve have 

g 

1 dp As 

p dT ~ BT 

Then using (2) in this relation ve find 


The molar heat of evaporation is 


As = R (1 + -r) 

v. 


L = TAs = HT (1 + gj) 

» § = \T\ ^ q » iT 


(3) 


1^ = 1+^ 


(f) From (3) m - * * M 

TaMng the logarithm of (2) , ve find 


+ 3_ = !n — = ln-S 


ha? 




where v is the volume per molecule of the vapor. 
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Hence 


-k~ = ln -£ 

^b 


(*) 


when v and v are evaluated at T, . 

gO D 


(b) In (4) we can take v to "be the volume per mole instead of the volume per molecule since ve 

I 

b 


have the ratio v /v . To estimate the order of magnitude of L/RT^ we let v^ = 22.4 liters/mole 

S ° -3 


and v q =» ^ 7ir o ^ 11^ = ^ it (.5 X 10 - ®) (6 x 10 2 ^) = 10 ^ liters/mole 


Thus 


L 


In - « 10 
10 _J 


i B 

V 

0 

V 

g 

mm mm 

vt, 

experimental 

water 

.018 liters/mole 

31 liters/mole 

.83 cal/gm°K 

1.45 cal/gm°K 

nitrogen 

.035 

6.4 

•37 

.64 

benzene 

.087 

29 

.16 

.27 


These results are remarkably good in view of the very simple theory considered here. 


CHAPTER 9 

Quantum Statistics of Ideal Gases 


9-1 


configuration 

no. 

of states 

0 

e 

3e 

MB 

BE 

FD 

XX 



1 

1 

— 


XX 


1 

1 

— 



XX 

1 

1 

-- 

X 

X 


2 

1 

1 

X 


X 

2 

1 

1 


X 

X 

2 

1 

1 


(a) 

00 

(<0 


= 

z be = 

^FD 


1 + e’ 2 ® + e" 6 ® + Se'® - 2,- 3 ® + Se' 1 *® 
1 + e‘ 2 ® + s' 6 ® + e‘® + e- 3 ® + 

e -® + ,-3»B «. .-"WB 


9.2 


(a) For an FD gas, we have 


-a-£E r _ x 

In Z = aN + Z In (l + e ) and n^ = — — 

r e r + 1 


Hence S = k [In Z + gE] = k 


-a-£E 

aN + Z In (l + e ~ ) + gi 
r 
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we find 


"but since N = Z n and E = Z 


_ _ n e , 

r r ’ 

r r r 


S = k 


_ -a-pe ~ 

Z n (a+ p€ )+ Z In (1 + e ' r ) 


r' r 
r r 


( 1 ) 


From n = — — it follows that 

r a+Pe 

e r + 1 


a + pe r = In (l-n r ) - In n r 

-a- P € r 5 r 


1-E 


( 2 ) 

(3) 


Substituting (2) and (3) in (l) we find 

S = k 


2 n r (lu(l-n r ) - la nj + Z la (1 + -H=-) 
r V sr r_ 

= -k Z Jn^ In n^ + (1-n^) In (l-n^)J 




(h) For the EE case 


-a-pe 


In Z = aN - Z lu (l-e r ) and n = — — 
v. r cc+Bi 


a+Pa 
e r -l 

The calculation of S proceeds in exactly the same way as in (a). We quote the result. 

S = -k Z n In n - (l+n ) In (l+n ) (5) 

r r r r r 

(c) From (^) and (5) we 6ee that in the classical limit where n^ « 1 

S -*• -k Z n In n 
r r r 

for both FD and BE statistics. 


9*3 


N / 

The partition function in £ /HI wher 


£ = Z exp 


K K K 

x y z 


6h“ r 2 2 2x 

- = (k + K +K ) 

2m v x y z ‘ 


( 1 ) 


The sum is over all positive k k k when we use the standing wave solution, equation (9.9.22). 

x y z 

The density of state is now greater than in the traveling wave case, and by (9.9*25) 

p d3 K= X-d3 K 

K ~ 7T 3 ~ 

If we assume that successive terms in £ differ little from each other, we can approximate (l) 
with integrals. 
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Thus 


„ _ 3^ 2 

Pm 

L e 


k =0 

X 


- 4 -^ L 

2m x / x , \ 

e ( — die ) 

' 7 r x' 


Hence 

which is identical to ( 9 » 10 . 7 )» 


3 

5 = (27T m IT ) 2 

■u-3 


9.k 


(a) The chemical potential can he found from the Helmholtz free energy, F = -kT In Z. 


-N 


F = -kT In |jy « -kT [N In £ - H In N + N] 


V / \2 

where we have used Sterling's approximation , In HI = H In H-H and £ = -=r (27r m kT) 


Then 




a) 


(b) The partition function for a single molecule of a 2 dimensional gas with binding energy -€ q 


00 


IS 


{;* = Z exp 


* * Be °iT exp 


t Ir^S^o 

* o f ~ | 

where A is the area. Thus £ ' = e £ where we put A = V 5 . 

II .*•< § 

Hence p' = -kT In £7 = -kT In — 


sir, 2^ 2* 

v — ( K +K ) 

2 m v x y ‘ 


N' 


Jar) 


( 2 ) 


„ dK dK 
2 x y 


(c) In equilibrium jj. = p. * and by (l) and ( 2 ) 


L.ill. 

N N' 


or 


V (2tt m kT)' 

N h3 


pe 


o A_ (27T m kT) 


Letting n* = N'/A and since p/lsT = H/V, we find 

-1/2 e As 

n' = ph ( 2 nm k^T' 3 ) e 


9 * 5 1 

(a) By equation ( 9 . 9 .IL) if L = L = L = L = we have 

x y z 
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( 1 ) 


/ 17 x 2rr 2 ii 2 '/ 3 , 2 2 2 \ 

e (V) = — - — (n + n + n ) 
r' ; m v x y z ' 


& . SL*! 4'% 

dV ‘ m \3 r x y z ' 


= — — ( 2 ) 

3 v 

(c) The mean pressure is p = Z n^ where is the mean number of particles in state r. 


— 2 y — 2 E 

P “ 3V r n r £ r “ 3 T~ 


(d.) The expression for p is valid for a system in which the energy of a quantum state is given 
by (l). For photons, / 1 


e = cfiK = ilc[2nV 3 [n 2 +n 2 +n 2 ) 
r V A x y z ' 


ne i e 

f rt 1 


„ / T\ J. T . — IE 

P r = _ 3 V 3114 p “ 3 V 


(e) The mean force on the wall is found by multiplying the change of momentum per collision, 
2mv, by the number of particles which collide with the wall per unit time, ^ nvA. Hence the 


pressure as 


Letting v 2 « v^ we have 


— i 2 

p = - (2mv) (g- nvA) = ^ rinrv ' _ 


- 1 2 2 N /I 2n 2 E 

TD = — nmv = — — ( — 221V ) = 

3 3 V *2 ' 3 V 


since N(= mv _ ) = S. Equation (3) is the exact result as shown in section 7.1^. 


(a) The dispersion is (Ap)^ = p“ - p 2 . By problem 9-5^ P = y — — g— , and it remains to find p 


o k 


P = 2 n p = — ^ £ n . e by (2) in 9-5 

r 9V r 


9 V 2 


Hence 


(Ap) 2 = (E 2 - 3 2 ) = (^E) 2 


— ^ 2 o 

(b) Since 2 = - rq In Z and (AE) = — ^ In Z, it follows that 

oS 
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(Ap) 2 = \ ^ In Z 

9v 83^ 


9V' 


2 53 


tut by (2) in 9-5, E = | V p. 


Thus 


(Ap)" 


2 Sp 
3V 53 


2kT 2 op 
3V 5T 


. _ v 

(c) Substituting p = — kT, we find 


-2 


(A P ) 2 /P 2 = (I M ) ^(Akt) 


N 

ST W 


2_ 

3N 


*9-7 

The probability that a state of energy e is occupied is proportional to the Boltzmann factor 
-e/kT 

e °. Since the vibrational levels have a spacing large compared to kT , we see that 
-e/kT 

e 0 is small and we may neglect this contribution to the specific heat. On the other hand, 

the spacing between rotational levels is small so that these states are appreciably excited. 

For a dumbbell molecule the rotational energy should be 

1 2 1 2 
e r = g A l“l + 2 A 2 w 2 

where A is the moment of inertia and w the angular frequency. We have neglected rotation about 
the axis of the molecule because the moment of inertia is small. At room temperature, we 

may use the equipartion theorem yielding = kT . Adding the energy of the translational 
motion gives for the total energy 


6 = e t + e r = 2 M o + M o = I M o 


Thus the molar specific heat is 


9-8 


0 ■ ”A It ■ I V ■ I R 

o 


The reaction proceeds according to 


and the law of mass action becomes 


2HD = H 2 + D 2 


\ n i>2 ^ h 2 \ 

- = K n = ~ 


“HD 




’HD 


( 1 ) 


From the results of section 9*12, we can immediately write down the quantities In the 
notation of that section. 
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V 2 P H 2 2A *2 

C = ^ [27r(2m) KP]^ e 2 2 

H 2 h 3 2¥l 


C = ^ [ar(2M) krr e 
2 h 3 


I *D„ ^ 


KT - | PfiWg 


I ps % 


5 m - \ ^ a' * ^ 

^ h 3 ii 2 

where we have assumed that molecules are predominantly in the lowest vibrational state and have 
treated the rotational degrees of freedom classically. 

The electronic ground state energies e 1 are approximately equal for the three species since 
this quantity is altered only slightly by the addition of a neutral particle to the neucleus 


(to form D) . Thus 


HD H„ 


S imi larly, the moments of inertia are of the form 

. _ 1 * _ 2 _ 1 ^2 2 
A 2 4 R o 2 nijHflig R o 

where R , the equilibrium distance between nuclei, depends primarily on electromagnetic forces. 
It is then approximately the same for the three species, and we have 

Aft a - 2 = = 1 (m^ ) 2 / 3) 

^2 ilg nD l 2 n 2'W ¥ mM 

Finally, the frequency for two masses coupled by a spring of constant k is 


ic(m^-m, 

m, m„ 
1 £ 


. Again, k must be the same for the three molecules, and 


k (m+M) \ 2 


HD V mM 


/2k\ 2 , 2m ^ 2 

' M ' ~ hn W “o 


/2k\ 2 /2M n 2 

111 ~~ HUM' U o 

Substituting the £ T s into (l) and using the results (2), (3) > and (2) yields 


K n * I $3 mT exp 


_l x 
m 2 +M 2 _ 

[2(mHM)] 


Since n^ = n^ , we have 
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2 


'N-: 


'“HD 


= K n*HD 


or 


^HD = 1 ^ 1 

N 1+2 JP 2 JF 




nM 


exp 


i I N 

iiio I . m 2 + MfL 

2kT 1 




9-9 

In a thermally isolated quasistatic process, there is no entropy change. Hence hy the first law, 

TdS = 0 = dE+pdV 

Substituting E = Vu and p = i u where u(t) is the energy density, we have 

0 = udY+Vdu+^udV 


r^k dV' 
'J v 3~ 


or 


du 

u 


T. 

X = — 

2 


r 1 dT * — k k 

/ k since u « T . Integrating, we find - — In 8 = 4 In — 

\J rp -*• 3 -*•-» 

i 


9.10 

(a) By the fundamental thermodynamic relation 

T dS = dE + p dV 

Then substituting E = Vu and p = ^ u, we have 

_ _ i _ k — — 

TdS=udV + V du + 1 u dV = — udV + Vdu 

But since u is only a function of T, du = ^ dT, and it follows that 


We also have 


dS = | | dV+^gdT 


dS = (If) dY + (t§) dT 

T V 


( 1 ) 

( 2 ) 


Equating coefficients in (l) and (2) yields 


/oS-v k u , ,oS\ V du 

o? j ' 3 T 8X1(1 dT ^ “ T dT 


(3) 


(b) Since (c 2 S/oV cl) = (o 2 s/gT o V), it follows from ( 3 ) that 
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du 4u 
dT ~ T 


or 


*/?=/# 


k - k 

Thus In T = In u + C where C is the constant of integration. Finally u « T • 

9.11 

The energy of the "black "body radiation in the dielectric is 

S = Vu = V 2- 

15 (c'h) 3 

where c' is the velocity of light in the material, c' is related to the velocity of light in 

vacuum by c 1 = c/n , n being the index of refraction. 

' o' o 


Thus 


E = V — ^1— n 3 = 2 

« (cH )3 o 


Here a is the Stefan Boltzmann constant, a = 7 r^k ^/60 c 2 fi 3 . 


Hence 


<r - <§ > - 


16 a n 3 VT 3 


( 1 ) 


Substituting v, the volume per mole, for V in (l) gives the heat capacity per mole, c^*. The 
classical lattice heat capacity per mole is cy = 3R. The ratio of these two quantities is 

c* 16 a n 3 vT 3 


C V 3RC 

For an order of magnitude calculation, we can let v = 10 cm 3 /mole and n = 1.5. At 300°K, we 
find Cy'/Cy ~ 10 ^ 3 . 

9.12 

(a) The total rate of electromagnetic radiation or power is 

p = (p(k rr r 2 ) = 47rr 2 oT' 

1l 

= (%)(10 2 )(5.7X10 _5 )(10 6 ) = 7.2X10 22 ergs/sec 

= 7.2 X lO 1 ^ watts 

(b) Since the total energy/sec crossing the area of a sphere of radius r* = 1km must be the 
same as that crossing the sphere of radius r = 10 cm, we have 

<p'(4 7T r' 2 ) = <?(% r 2 ) 

2 


or 


(?' =(?(fr) 2 = (fr-) 2 = (5-7 x io“ 5 )(id 6 )\^) 2 

r r 10^ 
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1 *] o il P 

= 5.7 X 10 ergs/sec cm = 5.7 X 10' vatts/cm 


(c) The power spectrum has a maximum at 


max 

kT 


he 


KT \ 


= 3 


max 


Thus 


9.13 


„ , x K) ' 27 )j? x VjP) . bj X 10- 6 cm - b 7 A 

(3) (1.4 x io" 16 )(io 6 ) 


(a) The power given off "by the sun is 

Jl O 

P = a T q (% R ) 

Half the earth's surface intercepts radiation from the sun at any time. Hence the rower that 
the earth absorbs is 


The pover that the earth emits is 


P = a T 4 (4ttR 2 )(^L) 
1 0 knL 2 


ll P 

P e = a T (km ) 


In equilibrium, the power incident upon the earth is equal to the power it radiates, or P^= P g 

k .. O. .On-* 2. k o. 

Thus 


oT k (knR 2 ) = tfT (% r 2 ) 

Wl 


T = T. 




R_ 

2L 


( 1 ) 


(b) From Eq. (l) 


9.14 


T = (5500) 


7 X IQ' 


.10 


-.1/2 


2(1.5*10 13 ) 


= 270°K 


We consider the liquid in equilibrium with its vapor (at vapor pressure p). By (7.I2.I3), the 
flux of molecules incident on the liquid is p/V2rr m KT' , and this must also be the flux 
of particles escaping from the liquid. Since jf cannot depend on the conditions outside the 
liquid, this is the escaping flux at all pressures . 

ThUS -V ~ /—■■■ " ■ ; - f r 

V27T m kT 

The mass of a water molecule is m = 18/6.02* IQ 23 = 2.99*10 23 . Thus for water at 25°C where 

) I p 

p = 23.8 mm Hg = 3*17*10 dynes/cm we have 


X 


3*17*10 


(2ttX2.99 X 10 - 23 X1. 38X10 _1 ^ 298)2 


= 1.14 x 10‘ 


.22 molecules 


sec -cm 
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9-15 

At the vapor pressure p(T) , the flux of molecules of the vapor incident on the wire is hy 
(7.11.13), <p = p/ *J2tt m hr'. When the metal is in equiUhrium with its vapor, the emitted flux 
must equal the incident flux, ® e = cp^. Then when the temperature is T, the number of particles 
emitted from the wire per unit length in time t is 


N = P(T) gnr t 

VarmS 

But N = £M/m where dM is the mass lost per unit length and m is the mass per molecule. 

Thus ^ * 

m V27T m hT 


or 


p(T) = ^ /§L 

; rt N 2nu 


where we have used u = mW^ and R = 


9.16 


v = 0 "by symmetry. 


By definition. 


J d 3 v f (v) 


J d^v f (v) 


( 1 ) 


At absolute zero, the particles fill all the lowest states so that the distribution f(v) is just 

1 2 2 ™ 

a constant and cancels in (l) . The maximum speed is at g mv~ = M* By symmetry, v 


1 2 
x 3 V 


and (l) becomes 


2 

v = 
X 


% 

3 


k a 
v dv 


krr 


2 , 
v dv 


_ 1 v 2 _ £ £ 
" 5 F 5 m 


9.17 

(a) 

At T = 0 all states are filled up to the energy |i. Hence 


the mean number of particles per state is just 1 and we have 

PM. 

E = / ep (e) de 
J o 


By (9.9.19) 


I i 


, oJL (2m) 2 

p(e)de = 2 p e 


1 2 x 3 
klT fi 


de 


where the factor 2 is introduced since electrons have two spin states. 

3 

v2 

\/ ( >m 

E 

J o 3m~ 


2 £ 
2 


v xaai. r %2 
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oo 2 f 

Since u = |L (3tt 2 |) by (9-15.10), it follows that 

E = | tyi (1) 

(c) Substituting the expression for p. in (l) we have 

2 2 

5 - Jj (37T 2 ) 3 £ (|) 3 N <*) 

As particles are added, each new particle can go only into an unfilled level of higher energy 
than the one "before. Hence the energy is not proportional to IT as it would "be if particles were 
added to the same level. 


9-18 


(a) Differentiating E in 9-17 (c) gives 


5 - - <i> 


2/3 £ a 5/3 

5 m V 


( 1 ) 


(b) Again from 9-17 (c) we find 

- _ £ I (3tt 2 ) ^ 
s ' 3 V 5 m V 

22 "5 

(c) For Cu metal, N/V = 8.1 x 10 electrons/cm 0 . Substituting into (l) where m is the mass of 
the electron, we find p ~ 10 5 atom. This enormous pressure shows that to confine a degenerate 
electron gas, a strong containing volume such as the lattice of a solid is required even at 0°K. 


9-1 9 


From the general relation TdS = ~i Q = CdT, we have that if the heat capacity is proportional to 

s> 

a power of the temperature, C « T 



T a_1 dT 


<x 


Thus the entropy is proportional to the same power of the temperature as is the heat capacity, 
(a) The heat capacity and entropy of conduction electrons in a metal is proportional to T. 
Hence the ratio of the final and initial entropies is 


ff _ iioo _ 

S. 200 

(b) The energy of the radiation field inside an enclosure is proportional to T^, and the heat 
capacity has a T dependence. 


Thus 


ff _ ,2000 

S i ~ '•1000 



7 1 



9.20 

(a) At T = 0, the total number of states within a sphere of radius k must equal the total 

F 

number of particles N. 

Thus 2 (| 7 r k 3 ) = N 

(2tt) 3 3 F 

1 

2 N ', 3 

“F = (3tt y ) 


or 


2 2 

Substituting p = fi k p /2m, we find 


fi 2 ,_2 N>, 3 

^ = ST (37r V } 


( 1 ) 


The number of particles per cm 3 , N/V, is N A p/p where is Avogadro’s number, p = .95 gm/cm 3 , 

and p = 23 . From (l) we find the Fermi temperature T_ = p/k = 36,000°K. 

F 

(b) To cool ICO cm 3 or 4.13 moles of Na from 1°K to . 3°K, an amount of heat Q must be removed. 

Q = vj c y dT 
2 

where is given by ( 9 . 16 . 23 ), = | R (y ^-) . 


Thus 


_ 2 , , p . 3 

Q = (4.14) R)(4t- — ) / T dT = -.0022 joules 

^ u ^ J 2. 


3 3 

Hence .0022/. 8 = .0027 cm of He J are evaporated. Note that a large amount of substance is 
cooled by a very small amount of helium. 

9.21 

2 A 

(a) If electrons obeyed MB statistics, the susceptibility would be x^ = Np m /kT where p m is 
the magnetic moment. In the actual case, most electrons cannot change their spin orientation 
when an external field is applied because the parallel states are already occupied. Only 
electrons near the top of the distribution can turn over and thus contribute to X. We expect 
then, that the correct expression for X is X^ with N replaced by the number of electrons in 
this region, 


= (?)N = (£-) N 
F 


eff ' p 'T. 


where T_ is the Fermi temperature. 
F 

Thus 


Np 


X = N 


m 


eff kT 


kT 


( 1 ) 


(b) Since Tp is on the order of 10 or 1C? , (l) yields a X per mole on the order of 10 
electrons obeyed MB statistics, X would be changed by the factor 


If 
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This is on the order of 100 at room temperatures. 


9.22 

When the field H is turned on, the energy levels 

of the oarallel spins shift ty -p. H while the 

m 

levels of the antiparallel spins shift "by +u m H. 
The electrons fill up these states to the Femi 
energy p as shown in the figure. The magnetic 
moment is determined "by the number of electrons 



Density of 
States 


which spill over from antiparallel states to 

parallel states to minimize the energy. This is (p m H)p(p Q ) where p(u Q ) is the density of states 
at the Fermi energy p Q . Consequently the magnetization is 


M = 2 ^ p(pj H/V 


But since p(p Q ) = ^ jj- , we have 




and 


y _ M _ 3 ^m_ 
B -2 „ 0 


, N 

where n = - 


9.23 

If we choose the zero of energy to he that of the electrons in the metal, the partition function 
of an electron in the gas outside the metal is ty ( 9 . 10. 7) 

£ - % (ar m Xlf e' VoAT 
ir 

where the factor 2 occurs because electrons have two spin states. The chemical potential p Q is 
ribund from the Helmholtz free energy. 

F = -hTlnZ = -kT [N In £ - NlmN + N] 

^o = 1 = 111 N = ^ ^3 M ) 2 6 VqAT 

The chemical potential inside the metal is just p = V Q - $, and in equilibrium the chemical 
potentials inside and outside are equal, i.e., p = p Q . 
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Thus 




flh- 


Hence the mean number of electrons per unit volume outside the metal is 

3 


n = ^ = i- (2?r m 1 st) 2 e - ' 5 ^ 


( 1 ) 


9.24 


Consider a situation in which the metal is in equilibrium with an electron gas. Then the emitted 
flux is related to the incident flux by 

* e = (1-r) «>. (1) 

The electron gas outside the metal may be treated as a classical ideal gas from which it follows 
that the incident flux is 


. 1-1 
= 5 nv = 5 


(2rm IT) 2 e“ $/<kT 


8 IT 
7im 


where we haved used (1) of problem 9.23. The emitted flux is then given by (l). But since the 
emitted flux cannot depend on the conditions outside the metal, this is also the flux where 
equilibrium does not prevail. Thus the current density is 


J = e® = v q - e ft~ r - ) (M) 2 


where e is the electronic charge. 


9.25 


There are (27 t) 3 d\ states per unit volume in the range (k : dx) and with a particular direction 
of spin orientation. Since the mean number of electrons with one value of k is + j] 


o 2 

where e = h k /2m, we obtain 


d 3 K 


f( K ) A i A- 


If we take the z direction to be the direction of emission, the number of electrons with wave 
vector in the range (jc:ds) and of either spin direction which strike unit area of the surface 
of the metal per unit time is 

4>(k) d 3 jc = 2f (jt) v z d 3 £ (l) 

Then the total current density emitted is given by summing (l) over all values of k which will 

allow the electrons to escape. That is, the part of the kinetic energy resulting from motion in 

2 2 

the z direction must be greater than V , or h /2m > V Q , where we use the notation of problem 
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9.23. In addition, ve assume that a fraction r of electrons of all energies are reflected. 
Thus the electron current density is 

r oo flkjfo) ** 


J=e* 

6 (&r) 3 




dK 


1 + exp §~( K X 2+K y 2+K z 2 ) “ft* 


where e is the electronic charge. The integral over is of the form 


xdx 


2 

-x -kz 


1+e 


2 

x -cz 


xe ’ = - i In [l+e“ x ■* a ] 


2 

-X -KZ., 

e +1 


h ““ j ■ » t ; fi 1 ** * exp 

Since V Q > p. and, in general, |g (h-V q ) I » 1, we can use In (l+x)= x to expand the logarithm 


and obtain 


T _ 2e (l-r) e 

fr ) 3 * 


P(n-V o ) „«2 2 m 2 2 


-pn k /2m n» -P» k /2m 


1/2 


dK / e 


dK. 


Each of these integrals is (2nm/pfi ) hy appendix A. b, and 


j = % me (l-r) ^^2 e -4>/bT 


9-26 


(a) Since the integrand in (9.17*9) is even, one can write 


p~ x 2 p« x 2 

-r / ex. „ / e x 

2 ” ^ -» (e X +l) 2 “ " Jo (e X +l) 2 


= 2 


e X +l 


+ 4 


- 1 0 


x dx 
o e X +l 


where we have integrated "by parts. Since the first term on the right vanishes 

— m 

x dx 


h- k l 


o e X +l 


( 1 ) 


(h) This expression can he evaluated hy series expansion. Since e <1 one can write 


-x 

x e x 


e X +l l+e" X 


-x -x, -2x -3x ! 

= e x [1 - e +e -e ... J 


= Z (-l) n e-^ n+1 ) X x dx 
n=0 


Hence (l) becomes 
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Ip = * S (-l) n f e “( n+1 ) x x dx = k 5 Jtil 


n 0 “ 


n=0 u o 

The integral is, "by appendix A. 4, eq.ua! to unity. Thus 


n=0 (n+l) J o 


e_y y ay 


I 4 z . 4 

2 n=0 (n+l) 2 


1 . ± + . 1_ 

C O p * * # 

2 2 3 4 2 


= 4s ( 2 ) 


(°) 










Fig. 1 Fig. 2 

One can use the contour C of Fig. (l) to write 


e 


Fig. 3 


dz = 2nl 


1 E - a 5 


Z Sin 77 z 


n=l n 7T 


n=l n 


or 


/ 


dz 


2 . 

z sin 77 z 


= -2i S 


(3) 


00 


Here the integrand -* 0 for z -»». The sum in (2) is unchanged if n -» -n. Hence one could equally 
well write the symmetric expression 

r dz 


J 2 ' 2l S 

g, z sm 77 z 

along the contour C’ of Fig. (l). Adding (3) and (4) one gets then 

T dz 


-4i S 


(5) 


(6) 


CHC 


, Z sin 77Z 


But here the contour can he completed along the infinite semi -circles as shown in Fig. Q since 

the integrand vanishes as |z|-»«. Except for z = 0, there is no other singularity in the 

enclosed area so that the contour can he shrunk down around an infinitesimal circle C surroun- 

o 

ding z = 0 as shown in Fig. (3) • For z « 0, 

1 1 


2 . 

z sm tt z z 


2 [7TZ - r TT^Z^ . . .] TTTT* 


1 r. 12 2 ■, 

= Ll - g77 z . . . j 


-1 


Hence 


1 r-, . I 2 2 

y [1 + g 77 Z . 

TVZT 


3-Z / ^ t \ /7T\ 


1 1 . 77 

' ] = “ + 51 + 

77Z 


2 . 

z sm 77z 


= (-2771) (g) = - 


. 2 
177 
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Thus it follows by ( 5 ) that 


3jj? = JUS 


or 


00 / ,\n : 

_ _ v (- 1 ) _ TL 

r\ 1 / 


n=l n 


12 


and by (2) , I = 7r 2 /3 . 


9.27 

(a) We expand e in the integral 


j(k) = 


,00 ikx . 
e dx 

- 0 = (e X +l)(e" x +l) 


x ikx , 
e e dx 

-= (e x +l) 2 


( 1 ) 


1 ® x , eo / \m 

e dx y (ik) m 

(e X +l) 2 m=0 mI X 


Interchanging the sum and integral yields 
\m 
m 


;<*> - l eg 

m=0 


X 00 

e m , v 
x dx = L 


(e X + l) 2 


m=0 




( 2 ) 


(b) Putting x=z in (l) to denote a complex variable, we evaluate the integral on the closed 
contour shown in the figure. Here C* is the path where z is increased by 27Ti and C" is a small 
circle around the singularity at i7T« Then the integral on the complete path is zero since there 
are no singularities in the region bounded by the contour. 



But the integral is 0 at x = + ® and cancels along the imaginary axis. Hence 


ikz . 
e dz 


ikz . 
e dz 


(e * +l)( e- 2 +1 ) 


C’+C" 


( e+ l) (e+i) 


( 3 ) 


where both integrals are evaluated frcm left to right. Along C' where z = x + 2iri, we have 


since e x+?ni = e X 


^ikz j 
e dz 


-f 


* g ik(x+2ni) 


1 ( e +l) (e +l) u -00 ( e +l) (e X +l) 


dx = e 


-27Tk 


ikx 


dx 


(e+l) (e +l) 
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The last integral is just J (k) . Thus from (3) we have 


(l _ e -2nk) f e 1 ^ 2 dz 


ikz , 
e dz 


(e z +l)(e" z +l) (e z +l)(e" z +l) 


(M 


C C' 

To evaluate the integral around the small circle C", we put z = 7ri + (; and expand the integrand 
in terms of £ 

r e ikz dz -V r 

J (e z +l)(e" z +l) 


C" 


_-7rh r e b _-nk f [l+ikC+...] 

e J ram dE ■ e J — — di 

c" c" 


= e 


-rck 


|cH-(27ri) (-ik)| = 27rke" 7rk (5) 

where we have used e = l+x+... and the residue theorem. Then combining (L) and (5) we find 

77 k 


^00 ikx . _ , -rrk 

e dx 27rk e 


^ -® (e X +l) (e~ x +l) l-e _27Tk 
(c) Expanding (6), we find 


sinh -n k 


(6) 


_ 1 _ W 2 + 

~ sinh 7T k - 1 + 350 


Then comparing coefficients with the series in (2) we can read off the I^'s. We find I = 77/3. 

9.28 

(a) According to equation (6. 9 .It-), the probability of state r of the system is 


P « e 
r 


-3E -cSS 
r r 


Then the mean number of particles in the state s is 

Z n g exp [-p(n 1 e 1 + n 2 e 2 + •••)- a(n 1 +n 2 +. . .) ] 


n = 


Z exp [ -3 (n 1 e 1 + n e + ...)- a(n 1 +n ? +. ..) ] 


1 2 


where the sum is over all n.,n , ... without restriction. 

Letting Z = Z exp [-3(11^ + n 2 e 2 +. . .) - a(n 1 +n 2 +. . .) ] (l) 


we have 


1 3 . 2 

n = - — ^ — In 

S 3 8€ s 


where the partial denotes taking the derivative where occurs explicitly. Evaluating the sum 


(l) in the usual way, we obtain 
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( 2 ) 


2 = ^E exp [-(cft^e^n^^E exp [-(a+pe^n^ 




Thus 


In £ = - E In (l-e 
r 


-a-Be 


and 


— 1 8 . 5 : 1 

" s = '^ £ s ^ a+j3e s 


- 1 


(h) Using the same argument which leads to (9*2.10) with the function 51 we obtain 


, 8n 
1 s 


(An r = - U s- 
v s' e Se 


where again the differentiation is performed where e g occurs explicitly. Thus 

a+Be 




- = E 2 (i + i-) 
(e S -l ) 2 


= n s (l+n s ) 

The correction term in (9*6.18) is absent showing that if only the mean number of particles is 
fixed, the dispersion is greater than if the total number is known exactly. 

(c) For FD statistics, the sums in (2) become 

r -(cf+Be-JI I -(ostpeJ] 

Z = [1 + e J L 1 + e J • * * 

~ -a-Qe 

In = E In (l + e r ) 


Thus 


and 


= i- 

3 0€_ ^ a+Se. 


+ 1 


(An) 2 = 


, on 
1 s 

P ci" 


a+3e . v 

e 5 - 2f 1_ > 

C£+3e 0 n s \ — j 

r s n \2 \n y 

(e +1) s 


= n (1-n ) 
s s' 
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9.29 


We can remove the restrictions on the sum over states by inserting the Kronecker 6 function 


Z = Z exp [-p(n 1 e 1 +n p e 2 +. ..)! B- 


2 N,Z n 


= Z exp [-?(n 1 e 1 +n 2 e p +.. .) ] f exp [ (K-Z n^) (a+ia' ) ] da’ 


where the sum is over n^ = 0,1,... for each r. Letting a = a+ia' and interchanging the sum and 


integral yields 


Z = fef 


£(a) da* 


where 


2(a) = Z exp [-(g+ge 1 )n 1 - (atf-£e 2 )n 2 . . . ] 


This sum is easily evaluated as in section 9.6 and 

In 2 (a) = -Z In (l-e r ) (2) 

r 

To approximate (l) we expand the logarithm of the integrand about a' = 0 where it contribute 


appreciably to the integral. 


In e^ 2 (a) = aN + In 2 (a) 


where 


8^ In 2 


= (a+ia')N + In 2 (a) + B^ia’) + ^ B 2 (ia') 2 + ... 
= ON + In£( a) + iCw+B^a’ - | B g a' 2 + ... 


2(S) =e a2I £(a) 


4B 2 a' 2 i(ir+B- L )a l 


As in section 6.8, we optimize the approximation by setting 


N + Bj =0 


h + = 0 


43„q: t2 


and (3) becomes e^ 1 2 (a) = e 0111 2 (a) e 2 

Substitution into (l) yields _12 


e °®£(a) r ~ i B 2 a ' 2 , , 
— e 

° -7T 


In z ~ aW + In 2 (a) = c fi* - Z In (1- 


-a-ge 


where we hs.ve used (2) . The value a is determined by (h) 
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IT + v- In a ) = IT - Z g— In (l-e r ) = 0 

oa t ' w Sa v ' 

r 


which, gives 


E 

r a+P€ r n 
e - 1 


CHAPTER 10 


fstems of Interacting Particles 


Since cr(u) du « and for spin waves m = Ak , it follows that 


cr(m) to « o 2 dm 


Then hy (10.1.20) 


P" JF™ v2 1 

“ / w (pbm) m 2 dm 


^ / 2 
J ° (e^-1) 


Substituting the dimensionless variable x = phm we obtain 


J72 _ x 


CO x 

e 


x 




Jo (e -1) 


The integral is just a constant not involving g. Hence 

cc p“3/2 k rp3/2 


(a) Ey ( 10 . 1 . 18 ) 


p OO 

In Z = PNt] -J In (l-e ^ lXD ) (^(m) dm 


where 


a D (°) = 


2 c J 

But since V = 6rr N ( — ) , vs find 


3V 2 

0-2 3 “ 

2tt c 


for m < 


for m > 


In Z = PNt} - / In (l-e S ‘ 5m ) m 2 dm 

“D 0 

In terms of the dimensionless variables x = p-nm and y = pbrn^, this gives 

In Z = ^ - ^ f In (l-e" x ) x 2 dx 


y 3 


y o 


This expression can be put in terms of the function L(y) by integration by parts. 
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^ 2 = y - pr 111 ^- e " x ) r " f 


r y 

o e X -l 


= y £“~ " 3H In (l-e y ) + HD(y) 

*"9 /t 

= gj- - 3N La (l-e ^ ) + ND (Sj/r) 


where k©_ = ttcu . 

D D 


(b) The mean energy is 


E = -Ni) + 


S.-^laZ.-fc^lnZ 


3”* <°p e r y 3j**y (_ £ 

l-e y y*” 'J o e X -l y^ \e y -:L 


= -Ntj + D(y) = -Nr) + 3H W D^yT) (2) 

(c) The entropy may be found from the relation S = k[ln Z + pE], Then from (l) and (2), 
have 

S = Nk [-3 In (l-e -y ) + lt-D(y)] 

- 0 n/ T 

= Nk [-3 In (l-e u ) + hD^^/T)] 


For y » 1, the upper limit in the function D(t]) may be replaced by «. 

Dw=%r# 

y^ J o e -1 
k . 

The integral is shown to be 7T /15 in tables. 

k 

Hence D(y) = — y y » 1 

5y J 

For y « 1, we can expand e ~ 1+x in the denominator 

D(y) = f = 1 y « i 

y 5 '-'o 

From the results of problem 10.2, we find for y » 1, or T « 0^ when we use y = = 0^/T 

l»z.g + 

58 d 3 
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k k 

= , T . 37r mr 
S = -Nr| + t — 

s =_^m(i-) 3 


For y « 1 (T » 0^) ve find, using e y « l + y 


In Z = § - 3N In 9 ^ + N 
3 = -Ifr] + 3H KT 

S = m [-3 In + 1+3 

10 . 1 + 

The pressure is given in terms of the partition function hy the relation p = ^ In Z where In- 
equation (l) of problem 10.2 


Thus 


where 


-0 /t 

In Z = g - 3H In (l-e ° ) + ND (e^ T) 


p , „ | . 3S g i ft— (- + ^2) + »KP 

l-e U 


aDtej/r)" 

"dfSj/r)' 

S(eJ/r)_ 

dV 

— — 


SDOSj/T) 


V T 


i^dx 3 

+ 


iD ^ s r/ T ^ (9 ^a7" 

5 - - ¥ *v« i 


e -1 
d0^ 


In terms of y = -(V/0. _) (d0.VdT) , this becomes 

k 3 p - » ^ + ^ 

For T « S D , D^T) -»y (|-) , and 


For T » 0 , DCSjy/T) -* 1 


p = N ia Sl; 7 . N KT 
^ 1 cV 5 "3 

CV ? V 0^ 




10.5 

Comparing the expression for the energy E, equation (2) in problem 10.2, with equation (l) in 
problem 10. 1+, we obtain 

5 . „ §a * , ^ 
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Thus 


(a) From electrostatics. 


/Sp\ k 7 /SE\ V 

“■ k( w> v -- ( S> t * K? T 


E « R 


(h) The induced dipole moment in terms of the polarizability a is 


P 2 = OE « E 


(c) The energy of a dipole in a field E is 


u = -P 2 *E = p 2 E « R 


By eauation (10.3.8) and (10-3-16), the log of the partition function is 

In z = In + N In V + | Ital (p) 

where n = W/V. The first two terms are the ideal gas partition functions while the last is the 
contribution of the interaction. Since the interaction is represented by an additive term in 


In Z, it must also be represented by an additive term in the entropy, S = k In Z - p 


le ^ z . 


Thus the entropy per mole is 


S(T,n) = S 0 (T,n) + S'(T,n) 


where S'(T,n) = k [§ n l(p) - P (§ n 1(B))] = n | l(p) - p ^ I (p)] 

and S'(T,n) = n A(t) - a [- § 3 2 ^ (^) 

To show that A(T) is negative, we must show that the derivative 0 / 0 B [l(p)/s] is positive. 

We have 

® ' hi" S’ « 

but -Bue" pu -e' Pu + l = 1 > 0 

e pU 

1 

since (l+Bu) is just the first two terms in the expansion of e p . Thus the integral is positive 


and A(T) is negative. 


A(T) = - I J jl -e _Pu ( 1+pu] hmR 2 dR 
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10.8 


From the arguments of section 10.6 we can write down the partition function for a two dimensional 


gas 


where 


N 


2 - i— (i!B) z 

- NJ < h 2 e > % 


o 

InZ^ = N In A + | ^ 1(3) 
1(3) = r (e' PU -l)27TR dR 


Denoting the film pressure by p , we have 




N 1 


N 2 


2 A 2 


or 


hT A 

W ’ “ + B 2 (T > 


where n = N/A and B 2 (t) = - ^ I. 


10.9 


(a) One can write down the mean energy from (10 • 7 • 3 ) and (10.7*5) 

S = -Ng|i H. S“ = -NS kDc B (x) 

^o i jz s' ' 

Here N is the total number of particles. The consistency condition is 

B sW = 2^S X 

We solve this equation in the given temperature limits using Fig. (10. 7.1) 
T « T 


(1) 


( 2 ) 


Here the slope, kT/2n JS, of the straight line is small and the intersection with B g (x) 

is at large x where B s (x) = 1. Thus (2) becomes 

, W 2n JS 

1 = ^ x or x = 


and 


2n JS 

a JE 
hT 


kT 


E = -NS 3jT (^i#) (1) = -2n NS 2 J 


(3) 


T = T 


In this region kT/2n JS is large and the intersection is at small x. We cannot, however, 
approximate B g (x) by the linear term ^(s+l)x (eq.10.7.12 ) since this would not cross the 
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line (l£r/2n JS)x. Hence ve must evaluate the second term in the expansion of B (x) . For 

s 

small x, ,3 

. . 1 x x . . , 

coth X-> — + -- T-T X « 1 

X 3 ^5 
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CHAPTER II 


Magnetism and Low Temperatures 


11.1 

The mean magnetic moment is 


M = r 


ZC-SEyBH) e 


-PE 


-£E 

Z e r 


1 S In Z 
P PH 


-PE r 

vhere Z = Z e 


11.2 


The change in entropy per -unit volume is given hy 

“ -/ ° <tl> ® 

J o T 


But from the fundamental thermodynamic relation dE = TdS - MdH, we have the Maxwell relation 


an 

c 

PT 

T H 


<t> - (?) 


vhere S and M q are the entropy and magnetic moment per unit volume. Since M q = XH and 


X = A/(T-s) ve obtain 


Thus 


SM 

, Ov AH 

^ST ' ~ ~ , \2 

H (T-0) 

AS = -f ° ^ dH = - — j -§- 
^o (T-e) 2 (T-s)~ 2 


11-3 

The change in the molar volume when the field is increased from H=0 to H=E is 

Av = r H ° (f|) dH 
J o T,p 


But from dE = TdS - p dV - McLH> ve see that 


© _=-© 

T,p T,H 


where v and M are the molar volume and magnetic moment. Since M = XH and X = a/[T-0 o (1-H2?) ] 


ve find 


S7 AS a H 

o 

^ T,H [T-S o (l-K*p) ]‘ 
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and 


Av 


H 

--L' 


o Ada Hflff 
o 


AS a 
o 


E 


o [T-S (Ha p)] 2 [T-e o (l-iarp)] 2 2 


11.4 

(a) From the relation dE' = TdS - pdV + HdM where E' = E - vfl 2 /8rr we obtain 

^V,T ' (ST V,H 

In the normal state, M = VXH ~ 0 since X = 0. Thus (gS/oH)^. ^ « 0. In the superconducting state, 

B = E + ^ = 0 


or 


— V 
M = - JE 


Thus M is independent of temperature and (SS/^iH) = 0 

v,T 

(b) Since the heat capacity is given by C = T(oS/cT) we have from (11.4.18) 


s _ S V ffl 
s n dT 


r -r - m fa. _c 1 . VTH /d^x 

C s C n 8T < S s S n' 5 tt 'dT' 47T ', 2' 

cur 


(i) 


where H is the critical field. 

(c) When T = T c , the critical field is 0. Thus 


c _c = 7Tc 
s n 4 tt 'dT' 


H.5 


From (l) in problem H.4 and E=E c [1 -(t/T c )~] we obtain 


VE 


3VE 


c -c = - 

s n 


2nT. 


T + 


27iT 


c jr- T 3 


But since C =aT + bT 3 , we have 


To satisfy the condition C s /T -> 0 as T ->0we must set the coefficient of T in (l) equal to 0, 
i.e., a = VE c 2 /2t r t/. 
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(a) The entropy is related to the heat capacity, dS = CdT/T. Then since at T=T c and T=0 the 
entropy of the normal and superconducting states is the same, we have 

m m 

a c C ffl a cCdT 
s(T = )-s ( „, _/ o ^r-/ o -V 


7 dT = 


From which it follows that 


a = 3 yh r 


(h) Since dE = CdT, the energy at T=T c and T=0 for the normal and superconducting states is 

T m 2 

p c 7T 

E n (T c) ‘ E n^ = ^ dT = -f- 

v o 


E s (T c ) - E s (0) = f Off 3 dT = 


But E n (T c ) = E £ (T c ) and by (l) a = 37A C • 


Eg (0) - E n (0) = - 


From the relation dS = CdT/T we have for the entropy in the normal and superconducting states 

T c 

S n (T c) " S n (0) = f 7dT = ^c 
J o 

T 3 

p C - ca* 

s s ( T c ) - S s ( 0 ) = J card's = — ~ 

o 

Since the entropy is the same for both states at T=0 and T=T , it follows that 0=37 /T . 


Hence at T=T 


% i 

c s " or 3 ‘ 3 

C 
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CHATTER 12 


12.1 

(a) 6 


Elementary Theory of Transport Processes 
(1) 6 (c) 6 

(*) l (c) £ 


12.2 

(a) t 

12.3 

(a) The equation of motion for an ion in a field £ is 

mx = e £ 


Hence 


By (12.1.10) , 


e£ t d , — e t z 

x = — T and x = — T 


t 2 ,J tV‘A 2E , 2 ,2 


and 


(h) The time to travel a distance x is 


- e£ 2 

X = T 

m 


t = -/§T ** = ^ T 

Thus the fraction of cases in which the ion travels a distance less than x is 

pvr 


i = 


T -t/r dt . -y/T 
e ' — = 1-e = . 


757 


12 . ^ 


The differential cross section is given "by 

a(n) = -2*. § 

where ca is the element of solid angle and s is 
the "inpact parameter" as shown. 

From the sketch, s = (a^ + a^) sin <p 





as 

d£ 


(a l + 


\ d® 

a 2 ) cos cp — 


Since ® = 90-9/2 and da = 27 T sin S' d9', it follows that 

ds _ "(a^g) sin g */2 
da % sin 9' 


93 



Hence 


(a +a ) sin 6 '/ 2 cos 9 '/ 2 

I - --- sva - ? 


■ i (v* 2 > : 


By arguemnts similar to those of section 12.1, one has that the probability that a molecule 


travels a distance x without collision is 


_ p- x /C- 


P = e 


where £ is the mean free path. Letting P = .9 and -t = ( V " IT na) yields 

-In .9 = x 1 / 2 ' na = x ^2' pa/kT 
Letting x =20 cm ana T = 300°K gives p = l.lt X 10 mm Hg. 


The coefficient of viscosity is 


l n v„t.^5-J2K (!) 

3 °o “ 

where we have used -t = ( n a ) To determine the cross section, we note that the volume 
per molecule is A/N^p where A = 39-9 gms, the atomic weight, and p = 1.65 gm cm , the density 
in the solid state. Thus the radius of the molecule is 


. m / 8 kT s 


„ ,3 _A_n 

^WT N, p 


and a Q = = 5*6 x 10 ^ cm 2 

From ( 1 ) we find rj = 1.1 X 19"^ gm cm 1 sec Then the ratio of the calculated value of tj to 
the observed value, 71 = 2.27 X 10 ^ gm cm ^ sec \is 


71 ,/tj = .U8 
'cal' ‘obs 


The coefficient of viscosity is proportional to T 2 since t] = mv/3 i/p' a Q ana v = (SkT/W) 2 . 

Thus if the temperature increas s, the terminal velocity decreases. The velocity is independent 
of pressure. 


The force on the inner cylinder is the pressure p times the area A = 2m RL. 

= pA = 2 ttRL 7] ^ 



We make the approximation that the air at the outer cylinder is moving with it, i.e., at 
velocity m(R+S) = cnR, and the air at the inner cylinder is stationary. 


Thus 


and the torque is 


Su cdR 
cir ~ 5 


G = ^R = IP . 

o 


(h) The viscosity is given "by rj = mv/3 VI" a . Since air is mainly composed of nitrogen, we 
let m = 28/6.02 X 10 2 ^, the mass of a nitrogen molecule. Putting ct q = lO -1 '’ cm 2 and T = 300°K, 
we find 

-4-1-1 

7] = 2 x 10 gm cm sec 


and 


G = 8 dyne cm 


a « vVc Z 
o 


12.9 

(a) We let 

where x,y, and z are to he determined. Letting L = length, T = time, and M = mass, the 
dimensional equation becomes 


But since C = FR S , it has dimensions MT 2 L S *^ 


Thus 


L 2 = L V X I^C 2 
t ‘ 2 L s+1 

L 2 = L x+(s+l) z t -x-2Zj^t+z 


which gives x+(s+l)z = 2, x+2z = 0, y + z = 0. 

On solving this system of equations, we find x = 4/l-s; hence the cross section a is proportional 
tovVf 1 - 5 ). 

{" 0 ) The viscosity is given hy t) = nrv/3c/2" o and here a ® 7^ 5 . Since velocity depends 


on 


the square root of the temperature, ve have 


^ c = ^+7)/2(s-l) 


12.10 

(a) Consider a disk of radius r and width Ax 

at the point x. There is a pressure difference 

2 

Ap across the disk giving force 7rr Ap and a 
viscous retarding force 27rr Axq ^ • Hence 

the condition that the disk moves without 
acceleration is 
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7ir 2 Ap = -2jtt Ax t] 


But since Ap = - Ax, where the minus sign is introduced because p(x) is a decreasing 


function of x, we have 


2 5p(x) 

7TT — §-* — u = 2ir TT\ 
OX 


Assuming that the fluid in contact with the walls is at rest, it follows "that 


du’ = — 


1 §p(x 


2t| ox 


, ^ 1 dp(x) / 2 

u < r ' x ) = 5^ -it 1 < r " a ) 


„ , . . , dp P 2 P 1 

For a liauid, w 2 - = — - — 
3 ox L 


The mass flowing per second across a ring of radius r* and width 


dr' is pu(r , )27rr'dr ' where p is the density. Then the total mass flow is 


M = / p^r^BHr'dr' (P^g) I ( a2 “ r ' P ) r,dr 


= 5 5 r (p i _p 2 ) 

(h) For an ideal gas, the density is p(x) = (p/RT)p(x) where p is the molecular weight. Thus 


the mass flow is 




(r' 2 -a 2 ) r'dr' 


= -Sfe p W 


Since M is independent of x, we obtain 

. rL 

M / dx = 


k 

7T pa 

5 t)BT 


which yields 


_ 71 ' a a /_ cL\ 

S tiRT L ^ P 1 " P 2 ' 


12.11 

Let T(z,t) be the temperature at time t and position z. We consider a slab of substance of 
thickness dz and area A. By conservation of energy, (the heat absorbed per unit time by the slab 
is equal to (the heat entering the slab per unit time through the surface at z) minus (the heat 
leaving per unit time through the surface at z+dz) . Since the heat absorbed per unit time is 
me dT/dt = pA dze dT/dt, where c is the specific heat per unit mass and p is the density, we have 
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by (12.4.9) k = — — v 

3^ ff o 

2 _s 2 2 

The specific heat per molecule is c = (3/2)h while a Q = 4 ttR = 4-7r(lO ) cm . The mean velocity 
is v = (83ST/7im)~ /2 . Choosing an intermediate temperature T = 286°K, we find v = 1.2 X l£p cm/sec. 
Substitution of these values into (l) yields 

Q = 1.2 watts/cm. 

(b) The thermal conductivity remains approximately constant until the pressure is such that 
the mean free path approaches the dimensions of the container, L = .6 cm. At lower pressures, 
the conductivity is roughly proportional to the number of molecules per unit volume n. Thus 
the heat influx is reduced by a factor of 10 when 

n = = lO 1 ^ molecules/cm^ 

10a/2 L c q 

From p = niff, we find p«3xl0mmHg. 
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12.14 

(a.) a 2 /a 1 = (^/rjg) (ug/i^) 1 / 2 

0-0 Kg/Vti — 

(c) D 2 /D 1 = t] 2 u 1 /ti 1 h 2 


(d) a = — ^— — r gives a = 1.1 X 10 ^ cm 2 , a = 2.8 X 10 ^ cm 2 . 

3^ 71 V ^ 1 2 

2 - 8-8 

(e) a = 7id“ gives = 1.8 X 20 cm, d 2 = 3*0 X 10 cm. 

12.15 

Tlie gases are uniformly mixed when the mean square of the displacement is on the order of the 
square of the dimensions of the container. We have, after N collisions, 

1 1 F J 


since the second term is zero hy statistical independence. 
Then 


,2 2 2 2 2 2 
C = V t=-VT 

* z 3 


2 2 2 2 2 2 
z = - v T n = 2 V ' rt 


where IT = t/c. Since t = -t/v = l/i/2 n a v, one has on neglecting the distinction between v‘ 


-2 


. 2 
and v , 


z 2 = J^U 

sM 05 


V = V| M foxi 

J j cp vir m 


o k 2 -1=; 2 

Substitution cf z =10 cm , a = 10 ' cm , and the given temperature and pressure yields 

1 

t as 10 sec. 

12.16 

A molecule undergoes a momentum change of 2(iV^ in a collision with the cube where jj. = raM/mtM, 
the reduced mass. Then the force is given by the momentum change per collision times the 
number of molecules which strike the cube per unit time, summed over all velocities. 


^ (H ) V2 * rv / f<vv): 


^2tt' 


d v 


We may approximate the exponent since v » V. Thus 
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? ( v .-V) 2 - ? V : 2 -ffmv- V 


(l+gnrv V) 


and 


l/2 p® _ §E v 2 
4. 0 /mBx 7 t 2 P 2 2 z 

Aot = 2tm ^ L J \ e 

w o 


(l+gmv^V) dv^ 


The first integral is just the force on a stationary wall as in section 7-l^* The resistive 
force is then 


l/2 _ p® 

= 2pn (g) L 2 / BmV v, 


gm 2 

- v 

3 _ 2 z 


27T' 


. I*n (g) lA I- L 2 

v 27r gm 


= pm v V L 

hy Appendix A. 4 and 7-H-13* Since we expect M » m, p » m, and the equation of motion "becomes 

dV — 2 

M ^ = -m v VL 


V = V exp 
o 


t i n v l2 *1 


The velocity is reduced to half its original value in time 

M In 2 


o = 


- t 2 
mnvl 


CHAPTER 13 

Transport Theory using the Relaxation-Time Approximation 


13.1 

By (13^.7) 


°el = ^ J 


Sg 2 

v 1 T V 

ce z 


2 2 2 

The value of the integral is unchanged if v or v is substituted for v . Hence since 

x y z 

2 2 2 2 
v = v + v + v ,we obtain 
x y z ’ 


el 


e" f ,3 3g 2 

— / d^V -r 2 - TV 

3 J - ae 


o 2 

The quantities t and cg/Se are functions only of v. Thus the substitution d J v = v dv sin 9 d9 dcp 
yields 


2 P 2ir p7T 


a , = — =— / dcp / sin 

el 3 Jq * J Q 


in SdS j' x v^ dv = - e 2 J" ^ r v^ dv 


(1) 
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13.2 

The function g becomes the Maxwell velocity distribution 

2 


5e - n ^2rr e 


3/2 - 5I_ 

23£T / 1 -v 


In terms of v = (PkT/m)^ 2 ,. this is 


$£ = -stt" 3 / 2 S£ 


"(V*) S 


*5 

v 


Substitution into (l) of problem 13*1 with s = v/v yields 


ne , \ 

a = (t) 

m ' 'a 


where 


13-3 


(t) = — / ds e“ S s' t(vs) 

3v7r ^ o 


Since e = \ mU 2 = in (U 2 + U 2 + TJ 2 ), we have 
^ c. x y z 


)U x cte oU^ cte U x 


= -m T d 3 U U 2 U t = -m 2 f d 3 TJ I s U 2 U 2 t 


TllUS T1 — -m I u. u u u * — — xu. / u. u v — u u 

1 J - a U x z X J - Se z X 

In spherical coordinates, the components of velocity became 


and 


U = U sin 0 cos cp U = U cos 0 
x z 


dij = I^dU sin 0 d0 dcp 


2 


^7^ 


o 2 P 5 

Therefore tj = -in" / cos“ ^ do I sin' 5 © cos^od© / U t 

Jo J o ^ o 


2 r W r . o1 . fl f” TT Sg tt 6 

= -m t T [cos 0 sin 0 - cos 0 sm 0jd0 / U t 


2 % P dg . 6 

= -m - j Q dU^J U t 


pi r p 27T 

But since / dcp / sin 0 d0 = it follows that 
J o o 


2 P 7T o 2?r 
dcp 

O J o 


r 1 = - / dcp j sin 0 d0 f dU || U 6 x = - 


m_ 

15 


( 1 ) 
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13A 


Substituting the expression 


|£ = V 2 ne " (U/ ^ - 

oe 


from problem 13.2 where v = (2hr/m) 1//2 into equation (l) of 13.3 yields 

ti = nJdT <t) 

r\ 

16 


where 


13- ! 




ds e s s^ t(vs). 


15 ^TT 


U/v 


Putting t — ( *72" n o q v) ^ into rj = hkT t and using v = (ShT/nm) 2 / 2 gives 

'fn “7 m KT 


n = T 


= . 1*3 


Vn ikT 


0 o 

This is larger by a factor of 1.18 than the mean free path calculation (12.3.18) and smaller 
than the result of the rigorous calculation (14.8.33) "by a factor of 1.20. 

13.6 


We assume the local eouilibrium distribution 


3/2 


/gnu r 1 D 2-, 

g = n(g^) exp - pmv ] 


( 1 ) 


where the temperature parameter 3 and the local density n are functions of z. Then in the 
notation of section 13 .3 > 

df(°) _ df^°) _ Sg Sn _ Sg S3 

Sn St ^3 St 


dt 1 


dt 


Sg Sn dz ^o) Sg SB dz ^c) 


-S 


^z dt 


}z dt 


Sg Sn 


Sg Ss 


v 


Sn 5z z z 

where the last follows since dz(t ) /dt = v . This excression is independent of f. Thus 

' o" o z 

( 13 - 3 * 5 ) becomes 


Af = 


df 


(o) 


dt' 


df 


or 


f = g - v 



+ |s 

S3 8 z 


( 2 ) 


Using (l) this expression becomes 
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1 Sn 3_ S§ mv 2 SB 

z n Sz 2B Sz 2 Sz 


It remains to determine Sn/Sz for which we use the condition 


v = - /d 3 v f v 
z n , / - z 


p ,3 f„3 2 1 Sn ^ 3 S3 mv SB 

or J d Z S v z -J d v v z Tg + ^ =0 

The first integral is zero since the integrand is an odd function of v . 


'3 2 

d v v z g 


1 Sn fz 1 _ mv 2 > |SS 
_n SS \2 B 2 )% z_ 


2 

where we have assumed t independent of v. By symmetry this equation also holds if v in the 

2 2 2 2 2 2 
integrand is replaced hy v x or v . Then since v = v x + v + v^. , we have 


fl Sn 3_ 

[n Sz 2g 

[l Sn 3_ 

[n Sz 2B 


'3 2 

d v v c 


m SB r ,3 t 
- w W— / d v v s = 


2 Sz J 


1, 2 
4 -2^v 


m SB 6 " I ^ 


dv = 0 


3 2 

on substituting (l) and d J v = v dv sin 9 d 9 do, and canceling common factors. The integrals 


are tabulated in appendix A.k 

fl Sn 3 _ Ml 1 r 2_v 5 / 2 m S3 15 /- ,2_> 1 

[n Sz 2B SzJ 8 ^ 3m' 2 Sz IS Bur 

1 Sn 1 S3 

or n Sz = 3 Sz 

Substitution of the equation of state n = p/kT = pB yields 

1 /'-n S3 „ Sp>, _ 1 S3 

(p SI + ? SF S Sz 

d P ^ 

or — = 0 

C.Z 

The pressure is independent of z. Then (3) and (r) give 

. r 2* 

. d3 5 2 iv 

f = s-sV SI |fe " 


1 ffl 5 1 mv 2 

T dz [2 " 2 hT_ 


Since there are no external forces and f is independent of time, the Boltzmann equation becomes 

Sf Sf f-f^ 0 ) 

I * Sv = V z Sz = " 
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where f ^ = g is the local Maxwell distribution of the previous problem. Letting f = g+f ^ 


where f^'« g we have 


z ciz 


8g Sn 


v - ^ = v + v 


[ 56 

z dp ^ 


.(1) 


Thus 


f = g-v T 


, 8g Sn . Sg S^. , 

L& ^ ¥ s£J 


This is identical to equation ( 2 ) of 13 . 6 . The rest of this calculation proceeds as in that 
problem. 

13.8 


The heat flux in the z -direction is given by 

Q„ = n <v (§ mv 2 )) = \ nm (v v 2 ) 


= i m / d^v f v v 2 
2.1 — z 


where 


f = g + 


gv T 


jz_ M 
T Sz 


5. 1 mv 

2 ~ 2 kT 


The integral over g is 0 since the integrand is an odd function of v . Consequently 

n r ’ 0 0 2 j 2 mv 


„ hit 1 5t 
Q z = T T ^ L n ^ 


3 „ „ 2 „ 2 2 ' 

2 . 2 


d v v v e 

— z 


{5-Bmv^j 


We can exploit the symmetry of the integral by replacing v^ by v / 3 - Then the only angular 

*3 2 2 

dependence occurs in d J v = v dv sin 9 dS dcp which yields wv dv when the integration over angles 


is performed. Thus we obtain 


o _ ZT 51 

^z " 3 T =z [fW 


3/2| r r- „ 

5 / dw e 

_ J o 


1 , 2 1«2 

6 " 2 ^ 8 • 2 ^ 

■ ~ -3m / aw e 

J o 


rhe integrals are tabulated in appendix A.V. 


TT 13T ST 
3 T Sz 


/iiO\ 


3/2| r 15 _ ,2 .T/2 


[ 5 -§ " %) ' -E» TT (fj) j 


9/2 


9 

5 hk T 


m 


ST 

Sz 


Hence 


5 nk 2 T 

k = o — IT - T 

2 m 


13.9 

The mean free path argument has 


1 — 0 k nk 2 T 

k = v ncv -c, = - T 

3 7T m 


since 


£ 


VT 


and v = (8kT/7im) 1//2 . Thus the constant relaxation time approximation yields 
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result greater than the mean free path argument by 


= I.96. 


13.10 

Prom problems 13-4 and 13-8 we find 

k = (3/2) (nk 2 T/m)T = 5 k = 5 £ 

T] n Iff t 2 m 3 h 

where c = (3/2)k, the specific heat per molecule. The simple mean free path argument gives 

il - £ 

r| — m 

Experiment shows that the ratio (k/t)) (c/m) lies between 1.3 and 2.5 (see discussion after 
(12.4.13)) so that the constant relaxation time approximation yields better agreement with 
experiment than does the mean free path argument. 


I 3 .H 

In the absence of a temperature gradient, the equilibrium distribution of a Fermi -Dirac gas is 


= e(c) ^ 1 


( 1 ) 


(27r) 3 fi 3 e P( £ “! 1 ) + 1 

In the given situation, f is independent of time but will depend on z because of the temperature 
gradient oT/ciz along the z direction. This temperature gradient alone would cause the electrons 
to drift, but since no current can flow, there is a redistribution so as to set up a field <£(z) 
to counter the effect of the temperature gradient and reduce the mean drift velocity v^ to zero. 
Thus the Boltzmann eauation becomes 


8f e £ (z) Sf 


f-f 


f z Si 


m 


5v 


We will assume that the temperature gradient is sufficiently small so that f can be written 
= g+f' - ', f' -/ « g where f 

3g , e£ og 


("0 (ll ( 2 ) 

t-f v ‘ ' « z where represents a small departure from local eauilibrium. 


,(l) 


Thus 


v 


z 3z m ov 


It is convenient to write this in terms of the dimensionless parameter 

x = p(e-ji) = P(w mv 2 -u) 


( 2 ) 


( 2 ) 


Thus 


8g 5g 5x 8g / \ 5 b x 5g 83 

5z _ 5x5z _ 5x' e ^'Sz _ BSx Sz 


8g _ 5g 5x 

' T 

Z Z 


sr = ^ sr = z If 
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Then from (2) we find 


f = g - v t !§■ 
z dx 


|1 dB „ e 

r r- X + 3 e fc 

'P dz H 


(4) 


It remains to find £ from the condition that v^ = i J d 3 v f -v^ = 0 where f is given hy (L). 

The integral over g is 0 since the integrand is odd. Letting t = t„ = const, for electrons near 

r 


the Fermi level, we have 




= 0 


or 


r a 3 v v 2 x + e 

P dz J — z dx H 


d 3 v v 2 
— z dx 


(5) 


The second integral can he expressed in terms of n, the total number of electrons per unit 

, _ . dg 1 dg 

volume. Since ^ = y— 

dx pmv dv 
z 2 


, 2 dg 1 

dv v tt 3 ' = — 
z z dx pm 


i . dg i r r 1 

' dv V nr 2 — = — V g] - — 

z z dv pm z° Bn, 

-00 z -oo 


g dv 


where we have integrated hy parts. Since g=0 at ± °=, it follows that 


'3 2 dg 1 r 3 n 

dvv ^=-— /dvg = - — 
— z dx 8m, — & pm 


(6) 


Proceeding to the evaluation of the first integral in (5), we note that x and dg/dx are functions 

2 2 2 *5 

only of v so that we can replace v hy v /3. The substitution d°v = v dv sin 9 d 9 dtp yields 


% v dv when the integration over angles is performed. 
Thus 


'3 2 dg % 

d~V v x v 2 - = -=- 
- z dx 3 . 


2 dg 
av v x ^ 


1 2 

Using x = p(— mv -p.) this becomes 


iH t±) 

3 'pm' 


5/2 


dx x 


(x4i3n) 3 / 2 


(T) 


From the properties of the Fermi distribution, we know that dg/dx = 0 except where e = 4 or 
x = 0. The integrand contributes appreciably to the integral only in this region so we can 
expand the factor x(x+p p) 3 ^ 2 about x = 0 to approximate the integral. Also, the lower limit -p p. 
may be replaced by -co with . negligible error. 

x(x+p^i) 3 / 2 « (Pi i ) 3 / 2 (x+ | 


and. 


2m 


3 


(2rr) 3 ft 3 ( e X +l) 2 


from (l) 
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Hence (7) "becomes 

ze { i .) /2 r 

3 j 


dx x(x+p 3 |_i) 


i/2 « 


2 3 / 2 u 3/2 n l/2 

2 ^ 

37 : P» 3 


(e X+ l) £ 


2 x 
x e 

-«(e X +l) 2 


These integrals are the ^ and. 1 ^ defined "by equation (9.16.9). It is shown there that 1=0 
2 

and I = 7T / 3 . Combining the results of (5), (6), and (8), we obtain 


(2pm) 1/2 1 Be _ ne £ 


3 Sz m 


and (^) becomes 


f = g-v t_ x 

z F dx 


n h 3 p 


Having found the distribution f, we can now calculate the heat flux 0 ^ 

Q z = n < v z (| mv ' 2 )^ = | 2111 < V 2 V ^ 

2 I. r ^ 2 

where (v v ) = — / d J v v v f 

' z n . / — z 


The integral over g is zero since the integrand is odd. 

Hence Q z = - ~f ± §f J d 3 v v^ 2 x|j - ± ■ ■ /d 3 v v/v 2 ^ (10) 

p 

The calculation of these integrals is similar to the steps leading to (8). Again v^ can be 

p 

replaced by v /3 and the integration over angles performed to yield The first integral in 


(10) becomes 


w r- - 1 - f «b> 


x(x+gn) 


5/2 


0 5/2 5/2 -1/2 r ® c 2 x 

2 ' n ' m ' / / S x ^ e 

” ‘ 3 H 2 B» 3 J - <X 2 ^ (e x +l) 2 

where we have expanded as before. The integral is I + ( 5 / 29 p)l 2 = 5 - n"‘ =: / 6 sn. 


Hence 


',3 2 2 Sg 5*2 3//2 4 3//2 m' 1//2 

d v v v x = - - — jr — - — — — - — 

S 9 g 2 ft 3 


The second integral in ( 10 ) is 


'3 2 2 Sg % f , 6 Sg 

d v v v -r 2 - = — / dvv r 2 
~ Z Sx 3 J D 8 x 
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u vyva 


3tt 


3ti- 


( 1 + 


5x 


2 ^(e X + l) 2 


dx 


Here we have neglected the squared term in the expansion since it is much less than 1. The 
integral is I q + ( 5/2 Pu) I. = 1. 

,3. , . m Obdi 


Thus 


d v v v r-^= 
z dz 


3 ?r 


ph- 


( 12 ) 


3/2 


O O — / 

Combining the results of (10), (ll), and (12) and using q = (fi /2m)(3rr n) , we find 


q = _ _Z i M 

^z 2 p dz 


5 2 JL_ + ^2 _n_ 

3 2 2 71 2_2 . 
^ m P m p J 


_ 7T 2 nT F 1 dp t£ nk^T_ dT 
^ "3 m T F dz 


3 m ^3 dz 


Thus 


7T 


2 nk. 2 ^ 


3 m 


13-12 

Combining the results of the previous problem and (13.4.12), we find 

"/‘.l = f ( i )2 T 

At 2T3°K this becomes K /° e j_ = 7-42 X 10 erg cm 2 deg ~ or in more standard units 
^ 1 

k/o ^ = 6.68 X 10 watt ohm deg . Experimental values at 273°K are 


element 

Ag 

Au 

Cu 

Pb 

Pt 

Sn 

W 

Zn 

K / CT el X 1C)6 
watt ohm deg 

6.31 

6.42 

6.09 

6.74 

6.85 

6.88 

8.30 

6.31 


14.1 


CHAPTER 14 

Wear-Exact Formulation of Transport Theory 


The mean rate of momentum loss per unit volume due to collisions is, by (14.4.22) 


(Ap) = / / / d 3 v d 3 v dfl* ff. Va (V,S’)( A ?) 


( 1 ) 


v v n 


where we assume central forces. We first perform the integration over angles by changing to 
center of mass coordinates. 

Thus (Ap) = m (v'-v) 
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"becomes • since v = c + *=• V and v* = c + — V* 

7 — — m— — — m 

(a e ) = n(v’-v) 


where p = j m “being the mass of the ions and m^, the mass of the heavy molecules. Since 

the mean values v and v are zero, it is clear that (An } = (Ap ) = 0. Hence we are left with 
x y 7 ' x y 

(Ap z ) = p(V z > - V z ) and 




dn' a(V,9 ')(Ap) 


( 2 ) 


V ’ is given in terns of V“ z and the scattering angles. Using the coordinate system illustrated 


in Fig. 14.8.1, we see that 


V ' = V'cos(V',z) = V cos (V',z) 


since V’ = V. Inspection of Fig. 14.8.1 yields 


V cos (V* ,z) = V 


cos(V, z ) cos S'+sin 9'cos cp' sin(V,z) 


= V cos 0'+ V sin S' cos cp' sin (V,z) 


and 


d(V 2 '-V z ) = pV z (cos 0'-l)+pV sin 9'cos cp' sin (V,z) 
2m 


o & i 

Substituting into (2) we obtain, since / cos cp’dcp’ = 0, 


>2n nir pm 

/ p(V '-V }o(V,0')sin 9'd9'd<j>'=2mpV / (cos 9' -l)a(v,0')sin 0'd9' 


o u o 


z z' 


For hard spheres, a is constant and 


p7T 

2mpV c / (cos 9'-l)sin 9'd9' = -47T pV a = -pV a,. 
z J 0 z z o 


where a,, is the total cross section. Equation (l) then simplifies to 


<AP Z > =p a t Jj V z V ff 1 d°v d-v x 


3 


Hi 


where 


3m 2 

■mS.3/2 - — 


(3) 


. 3/2 "3 S ( v 2 +v 2 +(v -u) 2 ) 

f = n(g) ' e 2 X y 2 

Here u is the mean z component of velocity of the ions. To avoid carrying repititious constants 


we let 


mg 


imjg 


cc = vf cs^ = — 


(M 
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f can be simplified by making the approximation that v » u, then 


„ 3/2 2 2auv 

/Q!\ ' _ -CCV z 


f = n0 e^ v e * n(^) ' e - ^ (l+2a u vj 


3/2 2 

x -av 


7T 


Equation (3) becomes 


- (Ap ) = nCTj.nn ( — ^ 


a^a n 3/2 nr -(ov 2 -kz,v 2 ) 

e (1+22 u v z )V z V d^v d^ 


( 5 ) 


77 

We again express v and v^ in relative coordinates 


1 

v = c + — — — V v. = c - — — V 

— — m+m^ — —l — — 


rhus 


2 2 x 2 2 (am -a m) (am^-tcnm 2 ) 

av + = (a-ta,)c + — — c-v + v 


m+m. 


On completing the square this yields 

(am -a,m) 




= (aw^) 


- + (af+a 1 ) (m+m^) - 


2 2 2 
(am 1 -a 1 m) 2 (am 1 + a-jn ) 2 

(af+ct^) (m+m 1 ) (m+m ) 


? aa-, 

= (a+a, ) iv + ■■ 

1 ' a+a^ 


where 


Also 


(am 1 - a^m) 

- “ - + (ct+a 1 ) (m+m^) - 


( 6 ) 


m (cm, - a m) m 

v = c + V = K - 7 — — va 4 — x V + — — V 

z z m+m^ z z (a-fa^) (rn-nn^) z m+m^ z 


a-, 


= K + 


z a-K 2 1 z 


(7) 


eve 


We wish to put the integral (5) in terns of V and the variable K defined by (6). It has been 
shown (14.2.6) that d^v^ d^v = d^£ d^V. In terms of K we ha^v 

d 3 c d\ = Jd 3 K d^V 

It is easy to check from (6) that J=l. Thus (^ becomes 

^ 2 pp -(a-ta 1 )K 2 ~a,+a V 

JJ 6 

The first two terns in the brackets give integrations over odd functions and hence do not con- 

1. C 

(on^) 


-(Ap z ) = nn 1 


a^a 


7T 


aa, 

|l+2a u K +2u — - . 

z a+a 1 zj 


V 


V V d 3 K d^V 

z — — 


tribute to the integral. Changing to spherical coordinates gives 

aa, „ 

1 2 
-~Z V c 

Q7+a n „5 


-(Ap ) = cr nn j ~~ ■ % 
z / t 1 (a-KLj 


5/2 r 

4?7 


-(a** ^k 2 

" iTdK 


^77 


'1 V^dV / 27T cos 3 sin 0 d9 
w o _ 


rhe integrals are found in appendix A. 4- and we obtain 
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/» x 2uu (° p j) 

p z ^ a t m l (a-tc^) 


5/2 r 


% ^ (a-ta^) “ 3// ^ 


/^L \ % 
LW / 3_ 


8 

3 ua t m ’ 


lm (i + i_)' 

1 m-Hn, N 7r 'm m . ' 




men 


where p. 


m+m. 


= - - ua x nn. 
8 1 


The momentum balance ecuation is then 


1 


which gives 


u = | 


ne£ + (Ap^) = 0 


£ 


Ca.ii V 2 u kT/m' 


Since j = neu = Cg^ £ , we have 


14.2 


3 -fn ne £ 


el 


T 


n i CT t V2|i it' 


(a) We will assume the distribution function 

f = f(°)(l-HJ>) 


where 


f (°) . nW (Sg M) 


3/2 


exn 


_ ^(z) v 2 


It is shown in problem 13-6 that the condition v^ = 0 yields the relation 

1 Sn _ 1 83 
n 5z p 

and that the pressure is independent of z. 

(b) The Boltzmann ecuation (14.7.12) becomes 

„(o) _ 8f ( °) , „ . 8f (o ^ 


+ v 


8f 


(o) 


= 


( 1 ) 

( 2 ) 

(3) 


at — dr m 8v 

Because the situation is time independent and there are no external forces, we are left with 

3 f (°) / */°) 


dz 


= V 


8n + 3 f^S3 _ f (°) 

z V n §z 2 3 8z 2 3z 


or (|g If) f (o) \ (l - an-fWf^oJvo a» (k) 

on using (2), (3), and (14.7.10). 

It remains to determine the function <J>. The left side of (4) leads us to expect that a good 
approximation will result if we choose 

* = Av (1-ay ) 
z 


no 



where A and a are constants. 0 mist satisfy the restrictions (it. 7-1*0 to he a good approximation. 

These become r , \ _ 

/ d^v f ' v (l-av~) = 0 

J ~ z 

fdW°) - v z ( 1- csv^) = o 

J ' d^v f^ v 2 v z (l-orv 2 ) = C 

The first and third integrals are zero because the integrands are odd. Similarly only the v^ 
component of v in the second relation will give a non-zero integral. We use this condition to 
determine a. That is. 



and 0 = A v z (l- ^ v 2 ) (5) 

We will use the procedure of section it . 8 to evaluate the constant A. Multiplying (t) by 4> and 
integrating over v yields 

k lt£ d3 X f<0) V& - f V 2 ) 2 . A fj" 6? Zl W f V 2 ) a[v 7 (1- f v z 2 ] 


The integral on the left is easily evaluated. We hav 

2 2Sm 2 


r * (o) 

/ d^v f ' ' 


= £r - + #) / dcp/ "cos 2 a sin 9d9 / n v° e 


- - /S 2 \ i • 

|V V V + ( — ) V V 

z ? z p Z 


( 5 ) 


(7) 


2 n27r pir 


pm 


'Sm 


■m3' 

2rr 


3/2 


3mv 


dv 


The integral of the second term in ( 7 ) vas evaluated previously. The third term is easily shown 
to he 7n/5f3ni. Thus the integral becomes 

I- (1-2 + J) = | S- (8) 

pm ' 5 ' 5 pm ' 

We split the right side of ( 6 ) into an integration over the angle and over v and v . We 
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define 


a\ a\ SE_) j 


X X x 


J = J ffl' V a a|j z (1- ^ V 2 ) 


Then using (8), equation (6) becomes 

(-l. M) (i 3-)- _EL_ M _ AI (o) 

^ Sz n 5 "2 5z 1 

171 

To evaluate J, we express A jv^l-av 2 )! in center of mass coordinates (we have reintroduced 
a = 3m/5). By (14.7.11), 

A Jv z (l-QV 2 ) =v z ’(l-av' 2 ) + V lz ' (l-OV 1 ,2 ) -V z (l-QV‘ £ ) - V lz (l-QV 1 2 ) 


Using v = c + i V > v^ = c - i V, we have 


(l-ov 2 ) + v. (l-av C )=(c + wV ) (l-ac^-ac'V - r V 2 )+(c - hr ) (l-ac“+ac-V - a r~) 

'j _LZ 1 Z X z x zxz 


= 2c - 2a c c 2 - ^ c V 2 - aV c»v 
z z 2 z z 


Then c = c* and |v| = { V 1 | gives 


A jv (l-av 2 )l = -a (V ' c-V' - V c«V) 
i Z z z 

Carrying out the product in (10) we find 


pcXT nTT (“ P ? 

J = -a do' sin S'de'cV c (V ’V ’-V V )+c (V ’V '-V V )+c (V ’ -V ) 

J o • J x v z x z x' y v z y z y' z v z z ' 

where the coordinate system is shown in Fig. 14.8.1. We notice that, except for constants, the 

first integral is (14.8.7). The second integral involves only a change of the index x to y hut 

is otherwise the same as the first. Consequently, by (14.8.20), these terms yield 


| a a W (c V + cV) 
2 Tj z v x x yy' 


where is given by (14.8.21). 
It remains to find 


r r ^ o p 

-a / dcp’ j sin S'dS’a Vc (V 1 -V ) 
On J n Z Z Z 


By (14.8.18) 


V,/ -V z ~ = V cos 3'+2V z V sin 9'cos 0'coscp' (|_-z)+ V sia 0'cos q'(f -z) 
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T dcp t (Y z ' 2 -v 2 ) = 2n V 2 (cos 2 0'-l) + ttV 2 sinV(t-z) 2 


To eliminate g, we write 


and (13) becomes 


z_ = 1 = (Syzj ‘ + (ij*z) + (£*z) 
*s2 , /O a.x2 

(l*z) = l-(l-z) 


dc?’(V ’ 2 -V 2 ) = 7TV 2 (2 cos 2 0 , -sin 2 9'-2)+rrV 2 sin 2 S’ 


The last term in J is then 


= (-37T V 2 + 7TV 2 ) sin 2 0' 


1 a o V cfv 2 - ~ 

2 T\ \Z 3 


and (ll) ana (it) give 


p 

|cmV(cVV + cVV + c V 2 - c rT) 
2 r| v xxz y y z zz z 3 ' 


We now evaluate 


= I “ V Cv a Cs-v) - c z 3 ) 

= r r a 3 v d 3 ^ f^) v^ci-ov 2 ) j 


Z It 


The evaluation of this integral is straightforward though tedious. In the center of mass system 


d 3 v d 3 v^ = d 3 c d^ by (14.2.6) 

,(0)^ (o) 2,m 6\ 3 pm , 2, 2\ 

; = n (— ) exp - (v +v ± ) 


2 ,mS\ J 

= n (— ) exp 

<c7i 


„ , 2 V \ 

-6m (c + -j7*j 


v„(l-av 2 ) = c^-csc^c 2 - ac„(c-v) - “ cj 2 + | V, - | V„ c 2 - | y_ c-V - g V 2 */ 


This expression is multiplied by J in the integrand 


„(o) „ (o) . . „ 

• i 7 1 _ 7 is clearly an even iimcn ion in 


c and V; hence we need keep only even terms in v (l-av^)j. This yields explicitly 

l.jpc ad r^f^d a cv) ^zSWSVl 

-c_ 2 I- + ?o. 2 c V + 55 e s V*+ I ■ : - Z ' r - 2 V 2 l (17) 


"z 3 3 Z 


Since f ^°^f 1 ^ has the simple four. (16) all these integrals are the standard ones encountered 
in this theory. We Guote the result of integrating over c. In (l3) the terms appear in the 
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X 


order that they appear in (17) 


I = 


2 a 


» 2 (f) 


31 


'■Hi 

2 ''pm' 


d^V 


SmV 2 


(18) 


W 2 - VV 2 - r V^V 2 - - V — + |2L y3 + 2_ y5 + « y3y 2 

z 2pm vv z ¥ v z 2 z 3 5jS v 12 £ vv z 


The integral hecomes on substituting a = pm/5, 

PmV 2 ,„ T 2 


d^V 

PmV 2 


T 

e a 


r W 


-^vV.vr + g v5 


V 5 a e 
I 


•>7T 


dV 2tt 


cos 9 7m6 ^2_ . _2. 1 


cos 0 - ^ + 


mPV 2 
60 , 


PmV 2 


dV 


- 1 <- 1 *> sy % 1 ~ T 

We let s = ^ V and 

= !(-« -e)/-=\% r 4s e_s2 


\Vgmj 


■ 2s \ 7 

a Hrr I s 


11 W, 


sin 0 d0 


102% n 
“ ' 15 (% 

vhere a is given by (12.8.26) . Substitution into (l8) yields 


T - 32 n~ 

25 7 ^F Z " 


and by ( 9 ) 


A = - 


- 4°) 


25 /HOT\ 1//2 8P 

5 TT ^ sj 

n _ 

i . _JS_ (y) l/2 v,(l- ¥ v 2 ) II 


32 no 


5 


8z 


To find the coefficient of thermal conductivity, ve evaluate the heat flux 

Q z = n < v z (| mv ' 2 )) 

= hn F d^v f^°^(l+<S>)v v 2 = T d^v f^°^$ v v 2 


1 . P , 3 -(o) 2 2,, mP 2 X 

= ^mA j d v f ' ; v^v(l — — v ) 


1 . 
= g 111 A 


_5 nS / 35 


(mp)' 
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Q - - 2 *! A « -= 5 - (JL) 1/2 M 

2 m 32 a V S 

7 


25 k_ hr kg 5g 
32 — Aj m 

n 

Therefore the coefficient of thermal conductivity is 


(c) By (it. 8. 32) 


K 


25 k I tt kg 
32 - a/ m 
7 


a 

1 


2 

= r- a 
3 o 


Hence 

it. 3 

(a) 


K 


75 k_ ) tt kg 
St a y m 


l7Tkg 

5 i/nmkg 

y m 

2t - 


a 

L TJ J 


-1 

15 k 
= TE 


3-75 


R 


(h) By (12. t. 2), the simple mean free path argument gives 

- = - = 1.50 - 

B m d 

Hence the Boltzmann equation gives a value for k/t\ greater than the mean free path calculation 
hy a factor of 2.5 


— X 10”^erggm' 5 aqg 1 

T) 

He 

A 

Xe 

Calculated 

1.55 

.780 

t 

OJ 

• 

Experimental 

1.55 

.785 

.250 


It.t 

Letting 



2-i 

pmv 

2 


Performing the integration over angles yields 



1= 2 
2^ r 


In n 




and by appendix A.t this is 
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Comparison with (9-10.10) shows that except for a constant, this is -S/k where S is the entropy 
per unit volume of an ideal gas. 

14-5 

(a) Differentiating the function 

H = J d 3 v f In f 

yields ^ = J"d 3 v || (in f +l) (l) 

Sf/5t can he found from the Boltzmann equation. We shall assume that there are no external 
forces so that the distribution function should depend only on the velocity and time. Then the 
Boltzmann equation becomes 

=//d 3 v 1 dfi* V0 (f^f - f x f) (2) 

Substitution of (2) into (l) yields 

ff = fff d3 Ii ® , Va(f 1 'f ‘ -f^ ) (In f+1) (3) 

This integral is invarient under the interchange of _v and v^ because a is invariant under this 
interchange. When the resulting expression is added to (3) we have 

f| = | fff d3 I d3 l! clfl' Va(f 1 'f’-f 1 f)(ln f^f+2) (4) 

Similarly, this integral is invarient if the initial and final velocities are interchanged because 
a is invarient under the inverse collision, i.e., 

a(v x v -» v^'v’) = ff(v^'v’ -»v 1 v) 

Thus || =| fff ^Zp' V'a (f^-f^f ')(ln ^'f+2) 

But V=V» and by (14.2.8), d 3 v‘ d 3 v^ = d 3 v d 3 v 1 

ft = if // d3 z d3 ii Va (v f i ,f,)(ljl f i ,f * + 2 ) (5) 

Adding (4) and (5) gives 

ft = ’ if ff d3 v d 3 v i dil* Va (f^f-f^Hln f^f'-ln f^) (6) 
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(b) Letting = y and = x, we see that since (in y - In x) (y-x) > 0, the integrand 

in (6) is greater than zero unless f^'f * = f f. Hence dH/dt < 0 where the eq uali ty holds when 
the integrand identically vanishes. 


Ik. 6 

The equilibrium condition is 


In f = Af B iv + B mv + B mv + C— mv' 
xx y y z z 2 


d 2 


(1) 


This can be written in the form 


In f = -A'(v-v ) + In C* 

o' 

where the constants A’, C’ and the constant vector v contain the constants in (1). 

-A ' (v-v ) 

Hence f = C' e 

v^ is clearly the mean velocity which can be set equal to zero without loss of generality. To 
determine A', we find the mean energy of a molecule. 

■“ 1 2 -A'v 2 


.3 1 2_ ,, / 2 1 2 

d v r iv f 47 r / v dVTir e 


- 2 


e — 


d 3 v f 


k r 2 a " A ' v2 

47T v dv e 


_ | V7T A r5//2 (| m) 


e = 


■R fZ/2 


3 m_ 
5 A’ 


Letting e = (3/2)kT, we find A’ = m/23£T. The constant C* is easily found from n = / d 3 v f which 
3/2 


gives C’ = n( ^a- ) showing that f is the Maxwell distribution. 
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CHAPTSE 15 

Irreversible Processes and Fluctuations 


15-1 

Let the friction force depend on a, v, and q in the manner 

f - aVV 

Then the dimensional equation where L = length, T = time, and M = mass becomes 

MLT' 2 - (LftHD- 1 ) V-V 1 )* - 

Thus x+y-z = 1, z=l, -y-z = -2. Solution of this system of equations yieldx x=y=z=l, and 

f = qav 
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15.2 

Since N k = R where is Avogadro's number and R is the gas constant, we have from (15.6.11), 

(x 2 > . JS- t - -2- t 


3TTTia 37niaR 


A 


or 


N 


RT 


A 37rna(x 2 ) 


Substituting R = 8.32 X 10 ergs mole ^ °K ^ and the given values yields = 7*0 X 10 2 ^ molecules 


mole . 


15.3 


(a) The energy of an ion in the field at point z is - eZ z. 


Hence 

(b) The flux is given by 


nfz+dz'l e Pe ^ C 2 ^ 2 ) 

^ - ~ peT - z ' - • 1 + ^* 2 


T „ 8n D[n(z+dz) - n(z)] 
J D " 8z * " dz 


(1) 


Then from (l) 

(c) 

(a) 

Thus 

15 .k 


= -nD pe£ 

j = nv = nu. £ 
d 


J + J = -nDp e / + mi£ =0 
D d 


H = £_ 

D M? 


S= F*(t) 


Substitution of v(t) = u(t)e 7 into (l) gives 

= ± e rt F*(t) 
dt m ' ' 


Hence 


u = u(o) + i / e 7 ° F'Ct^dt’ 


Since u(o) = v(o) it follows that 


= v Q e- Jrt f e 7t F*(t-)dt« 


( 1 ) 


( 2 ) 


15-5 

Ey dividing the time interval t into N successive intervals such that t = Nt, we can replace the 
integral in (2) of 15.I by the sum of integrals over the intervals of length t between t^ and 
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Then t' = kT + s where k is an integer and dt' = ds. 

N-l 

m 


Hence 


-v e'^ = e" 7NT ^ f T ^ kT+S ) F* (kT+s) ds 
° k=O m ^o 


= ^ e -7T(H-k) 1 r T e 7S pI(kT+s) ds 

k=0 m 


y S TT “1 

In the integral, s < t or e < e' = 1 since t is chosen so that t « 7 . It follows that 


v-v e" 7T = G, 

° k=0 k 


N-l 

= Z y k = y 

k=0 


where Gj^. = F'(kT+s) ds. 

15.7 

We integrate the Langevin equation 


^ = -yy + i F’(t) 
dt ' m x ' 


-1 

over a time t which is sma l l compared to 7 "but large compared to the correlation time t of 

the random force ^ 

v(t) - v(o) = -tv(o)t + i T F'CtOdt’ (1) 

^ o 

The time interval t is divided into N smaller intervals such that t = Nt. Then in the k 
subdivision t* = kT + s and dt' = ds. Hence (l) becomes 


v k = (1-7t) v k _ x + ijT F*[(k-l)T+s] ds 


We can find the velocity at t = Nt by an iterative procedure. Defining G k as in problem 15 .6 
have 

v X = fr-rO + G q 

v 2 = (1-7 t) v 1 + G 1 

= (1 -7 t) [( 1-7 t)v o + G q ] + G 1 
= (1-7 t) 2 v q + ( 1-7 t) G q + G 1 

Continuing this way we see that 

v N = (1 -tt) 1 * v q + ( 1-7 t) N_1 G o + ( 1-7 t) N_2 G 1 +...+ G^ 


Since 7T « 1, e = (e = (1~7 t) N it follows that 


N 


N . 


we 
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or 


\ ■ < S ‘ 7T '' -O + e ' rr(l " 1> 2 o * %+•••+ =„-! 

V N - V' 71 ' ** £ % 

* ¥ e- /T ^ N - k ) G . 
k=0 k 


since /t « 1. This result agrees with 15. 5. 
15-7 

G k =G 


= J ( F’(k+s)) 


ds = 0 


Hence 


and 


N-l N-l , \ 

Y - S y t . E 5 . 0 

k=0 k=0 

v - v e _rt = Y = 0 

o 


Then as t -»o, v -» 0. 

15.8 

From 15.5 ve have 


v = v e 
o 


■yt 


-w N-l — w N-l N-l 

Y 2 = Z y . 2 + Z S y. 


k=0 


* / j 


k y o 


It was shown in 15-7 that y k = 0. 

Thus ? . ¥ e- 2 "^) I 5 • 7 ¥ e 27Tk 

k=0 * k=0 

Where we have used Nt = t. The sum is just the geometric series and may he easily evaluated. 


e‘ 27t 

ri-e-1 

~ G 2 

l-e- 2rt 


l- e 2 7T 

~ Lr 

2rr 


2tt 2 

where e -1 = 27 T by series expansion since 7 T « 1. To find G we note that 


/ - 7 t, 2 ,.2 „2 l-e" 2rt 

(v-v o e . y . 0 27 T 


Letting t -»® yields, by equipartition. 


v 2 = M = G 2 
m 27 T 


Thus G 2 = ^ 7 T 
m 


and Y 2 = ^ (l-e” 2rfc ). 

m ' 1 
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Since (AT) 2 = Y 2 - Y 2 = Y 2 , it follows immediately from the central limit theorem that 

r ,11/2 f (v-v e'^) 2 ] 


T l1 l/2 

p ( v > t l v 0 ) = — ~ exp - 

Jbtl 2 ] 


2t&T( l-e"^) 


exp - 


, -7t%2 

m(v-v o e J 

2kT(l-e“ 2rfc ) 


15* ID 

From problem 15.5 we have 


G k = m/ F'Ckr+t’) at* 

J o 


,2 1 


m 2 ^ o 


at* / (f* (k*r+t *)f* (kT+t") ) at" 

^ o 


where we have changea the integration variable from s to t* to conform to the notation of the 
2 

text. G is inhepenhent of k allowing the choice of k=0 in (l) . 


g 2 = — 2 / at* / <F'(t*) f* ( t") ) at" 

m ^ o ^ o 


The integrana is the correlation function of F'vhich aepenhs only on the time aifference s = t"-t*. 
Making this change of variable, we have 

g 2 = — o r at* f (F*(t*)F’(t*+s)) as = ^ f at * F K(s)as 
J o 0 m'-'o w -t * 

where the last equality follows because (F'(t’)F'(t *+s)) = (F'(o)F'(s)) o = K(s). Since K aepenhs 

only on s, the integration over t' can be performea easily by interchanging the orher of inte- 
gration anh reaaing off the limits from the figure. 


g 2 = ^ I as / k(s) at 1 + / as / k(s) at 


* O p T 


1 T P T-S 


m I ^ -T w -S 


ds K(s)(t+s) + / ds K(s)(t-s) 


Since t » t* ana K(s) -» 0 for s » t*, s can be neglectea camparea 
to t ana the limits replacea by ± «. 



G 2 = ^ / K(s) as 
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By the results of 15 . 8 , G 2 = 213 ? rr/m. It follows that 

7 = afe r K(s) ds 


15.11 

(a) We integrate the Langevin equation 


f « -TV-ifft) 


over the interval t. 


r ( T ) 


dv = -7 F v dt' / F'(t')dt 


Since t is s mall , we can replace v hy the initial velocity v q in the first integral on the right 


and obtain 


Av = v-v^ = -7 V q t + i f F’(t')dt 


= -/V T + G 
o 


( 1 ) 


and 


Av = -7 v q t since G = C. 


(Av) 2 ' = (7 v q t ) 2 - 27 v q t G + G 2 


» G 2 = 2l3?7T/m 


( 2 ) 


where we have used the result of 15.8 and kept only terns to first order in t. 

(h) From (l), keeping terras to lowest order in t and noting that odd moments of G are 0 , we find 


(Av) 3 = -37 v t G 2 


(Av) = G 


2 lj. 2 3 4. 2 

Since G a t, G ' should he proportional to t . Hence (Av) and (Av) are proportional to r . 


O OO 

(c) From (2), (Av) = - 613? 7 t v /m. 


15.12 

From eq. (2) of 15 .h- 


(v(o)v(t)) = (v 2 (o)) e ^ e F'(t’)dt^ 

= (v 2 (o)) e _rt + ^ — (v(o)) f (e 7 ^ F'(t')) dt' 

'-I o 

But since F'(t') is random, the second term on the right is 0 . By equipartition 
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(v (o) } = 3£T/m 
Thus 


<v(o)v(t)> = f e-^ 


and since (v(o)v(-t)) = (v(o)v(t))., we have for all t 

(v(o)v(t)) = S e vW 


15-13 

Squaring equation (2) of 15 * ^ yields 

a -2rt nt . , 


2/ \ 2 -2/t e 

v ft] = v e + 2v 

v ' o o m 


e 7 ^ *F ' (t ')dt ' + e - ~~ f dt '[ e 7 ^ t ' +t ^F , (t , )F(t")dt" 


2 

m w o 


and v 2 (t) = v 2 e ~ 2yt + 6 - a f dt 'f e 7 ^ t ' +t ")(F'(t , )F , (t n )> dt' 


2 

in w o 


(1) 


since (F'(t')) = 0 . The correlation function in the integral is only a function of the time 
difference s = t"-t'. Making this change of variable yields 

f dt' f e 7(t’+t") (f ' ( t ')F ' (t ' ') >dt" = F dt ' /^ ds e 7 ^ 2t,+S )K(s) 

J o j o O ^ -t ’ 

where K(s) = (F'(t ')f '(t '+s) ) = (F '(o)F ' (s) ) . The evaluation of this integral is similar to 
that of problem 15 . 10 . Thus interchanging the order of integration and reading the limits off 


the diagram of problem 15.IO with t substituted for t yields 
\t—t' _ nO ot 


.-7 


I dt'/ ds e 7 ^ 2t +S ^K(s) = f ds / dt 'e 
o J -t ' 


-t d -s 
> o 

-t 


( 2 t ,+s )j^( s ) + I" ds T dt 'e 7 ^ 2 ^ ,+S ^K(s) 
O J o 


vs » r e 27t - e " 27S ' 1 r t ys r e 2 Ht-s)_iS 

ds e /S K(s) 1 M- / ds e 7 K(s) 1 <- 


27 


Since K(s) is only appreciable when 7s « 1, we can put e 7S w 1 in the above and find 


ds K(s) 


(° gn -l) 


2 7 


r\ 

+ / as 


K(s) hfgni = f «.*.)-*• iilil 


j 2rt -i) r t 

2 7 J. 


,27t 

27 


K(s) 


Substitution in (l) yields 


! (t) = v 2 e' 2rt + f ds K(s) 


v 


27 m 


as t-s-co^ v («) = kr/m by eauipartition and 


v (-) = ^ = -A- / ds K (s) 


Hence 


7 = 


Pm VP 


27m 


ds K(s) as in (15.8.8) 
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15.1^ 

(a) We will denote the Fourier coefficient of F'(t) "by F'(cjd) and that of v (t) by v (m). That 


F'(t) = f e^ 0 F'(m) dm and v(t) = f e^ 1 ^ v(w) du 


Substitution in the Langevin equation gives 


iwv(m) = -yv(as) + — F'(m) 


or v <®) = m( 7 iiu) 

Denoting the spectral densities of F' and v hy J (as) and J (os), we have from (i 5 .i 5 .ll) 

x V 

J f ,(cd) = | |F’(m)| 2 =| |m(r+4fl)) v(cd)| 2 


= m 2 ( 7 2 -to 2 ) ^ | v (as) | 2 


(b) Ey (15.10.9) 


J p .(®) = m 2 (-/ 2 -to 2 ) J v (co) 


K(s) = (v(o)v(s) ) = ^ e" 7,si 


J v (co) = |p f K(s) e' 1 ^ ds 


; (-ixD+ 7 )s ds + r a e (-ios- 7 )s. 


T / \ _ KT 1 1 , 

v' ' 27 im (7-iu) (-7-im) 


KT 27 


2nm 2 2 

L 7 --to . 


(c) Substitution of (b) into (2) gives 


«*■=.(“) = J mKT 7 


or J f / + ^(cd) = 2J F (m) = | mfcT 7 

This result depends on co « (t*) where t* is the correlation time of the fluctuation force, 
since in deriving (3), one used a time scale such that t » T* . Consequently, cdt =1 » un*. 



